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Abstract 


An experimental investigation of the effects of angle of attack on hypersonic 
boundary-layer stability on a flared-cone model was conducted in the low-disturbance 
Mach-6 Nozzle-Test-Chamber Facility at NASA Langley Research Center. This unique 
facility provided a “quiet” flow test environment which is well suited for stability 
experiments because the low levels of freestream “noise” minimize artificial stimulation of 
flow-disturbance growth. Surface pressure and temperature measurements documented the 
adverse-pressure gradient and transition-onset location. Hot-wire anemometry diagnostics 
were applied to identify the instability mechanisms which lead to transition. In addition, 
the mean flow over the flared-cone geometry was modeled by laminar Navier-Stokes 
computations. 

Results show that the boundary layer becomes more stable on the windward ray 
and less stable on the leeward ray relative to the zero-degree angle-of-attack case. The 
second-mode instability dominates the transition process at a zero-degree angle of attack, 
however, on the windward ray at an angle of attack this mode was completely stabilized. 
The less-dominant first-mode instability was slightly destabilized on the windward ray. 
Non-linear mechanisms such as saturation and harmonic generation are identified from the 
flow-disturbance bispectra. 
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1 Introduction 


The development of efficient hypersonic-flight vehicles of the future will require the 
accurate determination of the aerodynamic forces and heat-transfer rates which are 
encountered. Examples of such vehicles include atmospheric re-entry spacecraft, 
aerobraking orbital-transfer vehicles, high-altitude aerospace planes, and intercontinental 
ballistic missiles. The performance and design of these vehicles strongly depend on the 
state of the boundary layer. The boundary layer is the region of flow very close to the 
surface where viscous effects dominate the flow physics. This viscous-flow region 
naturally evolves from an orderly laminar state to a chaotic turbulent state. The local 
aerodynamic-friction force and heat-transfer rate increase considerably when this laminar- 
to-turbulent transition occurs. Furthermore, peak-heating rates are often attained within the 
transitional region. Thus, predicting the region over which these higher aerodynamic and 
thermal loads occur is of crucial importance in determining the overall vehicle loads to 
accurately evaluate the vehicle performance. 1 

The motivation for conducting boundary-layer stability experiments is derived from 
the need to validate numerical prediction methods. Only through validation with reliable 
data can confidence be placed in numerical results. Also, stability experiments seek to 
increase the understanding of flow physics. Ultimately this new understanding may lead to 
improved mathematical models which will extend the capabilities of the predictive tools 
used in hypersonic vehicle design. Angle-of-attack effects on a conical model are studied 
for the following practical reasons. The conical geometry represents a general 
axisymmetric forebody or a slender vehicle. In addition, angle of attack results in transition 
asymmetry which may significantly influence the dynamic stability of hypersonic vehicles. 2 
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1.1 Background 


In a low-level disturbance environment, and when bypass mechanisms are not 
present, the laminar boundary layer will naturally evolve into a turbulent boundary layer. 3 
The natural-transition process is described by the following sequence of events. 4 First, 
both body and ambient disturbance fields are generated. These initial disturbances are then 
modified by the body flow field and are subsequently internalized by the body boundary 
layer. This internalized disturbance field experiences linear amplification within the laminar 
boundary layer via normal wave-like modes. When these disturbances grow sufficiently, 
non-linear interactions and spectral-broadening effects take place. Finally, the onset of 
transition is observed by the first appearance of localized turbulent regions in the flow. In a 
high-level disturbance environment, however, the above process may be bypassed, and 
transition would occur through other processes. 

1.1.1 Theoretical & Computational Tools 

The most significant portion of the natural-transition process is the linear- 
amplification stage of normal-mode disturbances. Linear Stability Theory (LST) is a useful 
tool for predicting such amplification at a given streamwise location within a specified 
laminar boundary layer. LST is valid within the parallel-flow approximation, which 
assumes small mean-flow variations along the streamwise and span wise directions. A 
normal-mode disturbance of the form 

<D(y) ei a(x-ct) (1.1) 

is superimposed on the laminar mean-flow solution, where represents a physical quantity 
such as velocity, pressure, temperature, or density. This wave-like disturbance propagates 
along the streamwise direction, x, and y is the distance from the surface. The boundary 
layer is specified by the mean-flow profiles along the y direction of the velocity, pressure, 
temperature, and density. The flow disturbances of primary interest here can be classified 
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as either first or second mode. First-mode disturbances are typically low-frequency 
velocity fluctuations. These vortical disturbances dominate the low-speed transition 
process. In hypersonic flow, however, second-mode disturbances are also present and 
may in fact domi nate the tr ansi tion process.** These second-mode disturbances are typically 
high-frequency fluctuations in pressure. 

For the special case where the phase velocity, c, is real and the wave number, a, is 
complex, the latter can be written as 

a = 0Cr + iai. (1-2) 

Substituting Eq. 1.2 into Eq. 1.1 above, the imaginary component of a produces a real 
coefficient of amplitude A = e-cti(x-ct). This real coefficient provides the mechanism for 
linear normal-mode amplification. When the term (-otj) is either negative, zero, or positive 
the disturbance is either stable, neutrally stable, or unstable, respectively. A local stability 
analysis can be applied at a series of streamwise stations along the body to evaluate - 0 C|(x). 
For parallel flows, the integrated-growth factor, N(x), can be evaluated by 6 

N(x)-N(xo)=| -a,-© = In (1.3) 

Jxq 

where xq is a reference location, typically chosen as the upstream point of neutral stability 

where -ai = 0 such that N(xo) = 0, and £ is the variable of integration along the streamwise 

direction. The above expression also shows the relationship between the growth factor and 

the disturbance amplitude. The local rate of change in amplitude is given by 

— = Ax-oti (1*4) 

dx 

thus, -ctj is commonly referred to as the spatial-amplification rate. 

A method for LST-based computation of amplification rates and integrated-growth 
factors for the normal-disturbance modes is described by Malik.? The application of this 
method, termed the e N method, to transition prediction requires correlation of integrated- 
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growth factors with transition measurements. To successfully apply the e^ method, 
however, a reliable database of transition data where body geometries and flow conditions 
are similar is needed. 4 This necessity provided the motivation for previous transition 
experiments on cones at angle of attack in hypersonic flows. 

1.1.2 Transition Experiments 

Several experiments have documented angle-of-attack effects on transition for sharp 
cones in hypersonic-flow conditions. 8 " 14 These experiments consistently observed 
transition shift downstream on the windward meridian and upstream on the leeward 
meridian, relative to the zero-degree, angle-of-attack case. A few of these studies have also 
documented the circumferential transition patterns over regions where the cross flow 
velocity is significant 8 " 10 * 12 Only the location of transition was documented in these 
experiments. Stability experiments, on the other hand, investigate details of the instability 
mechanisms which lead to transition. 

1.1.3 Stability Experiments 

A comprehensive review of earlier hypersonic boundary-layer stability experiments 
was given by Stetson and KimmeL 15 This report described the application of hot-wire 
anemometry diagnostics to measure amplification and growth of the normal-mode 
disturbances. Previous stability measurements 1 ^ and LST-based computations 1 ^ have 
confirmed the dominance of the second-mode disturbance for a sharp cone at a zero-degree 
angle of attack. A significant discrepancy between the measurements and computations 
was noted in the peak-amplification rates, however. 15 These variations may be attributed 
to the influence of the freestream disturbance field in the test facility, which is not 
considered in LST-based computations. Furthermore, the greatest discrepancies between 
the stability measurements and LST predictions corresponded to those regions where non- 
linear disturbance interactions were strongest 18 The measured growth of high-frequency 
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disturbances within the range of stable frequencies predicted by LST was attributed to 
harmonic generation of the dominant, second-mode instability. Additional non-linear 
interactions between freestream disturbances and dominant-instability mechanisms may 
occur, however, where the freestream-disturbance levels are high. 

Stability measurements on a sharp cone at angle of attack in a hypersonic flow were 
first presented by Stetson et al. 19 Significant findings of this study are summarized in the 
following. The amplification rates of second-mode disturbances were not greatly affected 
by angle of attack. The region of amplification shifted, however, downstream for the 
windward ray and upstream for the leeward ray. The second-mode frequency increased on 
the windward ray, scaling with the thinner boundary-layer thickness. On the leeward ray, 
the second-mode frequency decreased, obscuring the first-mode disturbances. Growth of 
the first mode was reduced on the windward ray. These angle-of-attack experiments were 
conducted in the same wind-tunnel facility as the sharp-cone experiments at a zero-degree 
angle of attack.! 6 Freestream disturbances in such a test facility may artificially stimulate 
growth of the normal-mode disturbances and prematurely induce transition as discussed 
below. 

1.1.4 Wind-Tunnel Disturbances 

The primary disturbances in a supersonic or hypersonic wind-tunnel facility 
illustrated in Fig. 1.1 are of three types: entropy, sound, and vorticity modes. The 
entropy-mode disturbances arise from temperature spottiness and are negligible if the flow 
is thoroughly mixed prior to entering the test section. The sound mode propagated from 
the valves upstream of the settling chamber are eliminated through the use of fine meshes in 
the settling chamber. These mesh screens also reduce the amplitude of the vorticity-mode 
disturbances. Remaining vorticity disturbances are largely attenuated by the expansion 
downstream of the nozzle throat. Most importantly, however, eddy structures in the 
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turbulent boundary layer on the nozzle wall radiate acoustic noised These disturbances 
propagate along Mach lines from the nozzle wall towards the centerline axis. These 
acoustic disturbances, or pressure fluctuations, can stimulate disturbance growth in the 
boundary layer of the test model which may prematurely induce transition to turbulence, as 
illustrated in the figure. 

The freestream disturbance field must be considered when comparing transition 
measurements from different test facilities and extrapolating these measurements to 
atmospheric-flight conditions. 21 The above collection of transition experiments 8-14 were 
conducted under a wide range of flow conditions in a variety of test facilities which include 
a b allist ics range, 14 a shock tunnel, 13 and wind tunnels of different geometries. 8-12 
Variations in the transition measurements are thus expected due to the unique disturbance 
environment of each test facility. In addition, hot-wire measurements obtained in stability 
experiments include contributions from both the normal-mode disturbances in the model 
boundary layer and the facility disturbances. Small-amplitude normal modes may be 
obscured in a high-level disturbance environment Therefore, to increase the signal-to- 
noise ratio of the normal modes for stability measurements, the facility disturbances should 
be mi nimize d The above deleterious effects of facility noise on stability and transition 
measurements prompted the development of low-disturbance wind tunnels, described in the 
following. 

1.1.5 Low-Disturbance Wind Tunnels 

To provide a suitable environment for stability and transition experiments, NASA 
Langley has developed a series of low-disturbance wind tunnels. 22 Significant features of 
such test facilities are illustrated in Fig. 1.2 and discussed below. First, the air supply is 
highly filtered to remove particles from the flow. A series of progressively finer mesh 
screens significantly reduce the vorticity and entropy disturbances in the settling chamber 
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producing a very uniform flow. The dominant acoustic disturbances are reduced by 
delaying boundary-layer transition on the nozzle wall. This is achieved by applying the 
following techniques. Wall suction is applied through an annular bleed slot upstream of the 
nozzle throat to remove the turbulent boundary layer on the settling chamber wall. A new 
laminar boundary layer is established on the nozzle wall at the bleed-slot lip. The nozzle 
surface itself is highly polished to delay surface-roughness induced transition. In addition, 
a straight contour is used just downstream of the nozzle throat to delay the development of 
Gortler vortices. Because transition on the nozzle wall is delayed, the level of acoustic 
noise radiated onto the test model is reduced as illustrated in the figure. The present 
stability experiments were conducted in such a low-disturbance facility. The specific 
objectives of the present study are outlined below. 

1.2 Objectives 

The primary objective of the present work was to conduct the first hypersonic 
boundary-layer stability experiments on a conical model at angle of attack in a low- 
disturbance wind tunnel. The stability measurements obtained in the low-disturbance 
environment are suited for validation of stability theory prediction methods. The present 
stability measurements are also unique because new hot-wire anemometry diagnostics were 
applied. An additional objective was to generate laminar-flow computations for benchmark 
comparisons with mean-flow measurements. These mean-flow measurements were 
obtained in support of the stability measurements. The following chapter details the 
approach taken in the present experimental and computational work. Following this, the 
results are presented and discussed. Finally, this work is concluded with an overview of 
the significant findings and recommendations for future study. 
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2 Approach 


The experimental apparatus and procedure implemented in the present study are 
described in this chapter. This experimental work was complemented by a computational 
study of the laminar mean-flow boundary layer on the flaied-cone model at angle of attack. 
The computational method employed is also described in this chapter. Finally, post- 
processing techniques applied to the experimental and computational data are described. 

2.1 Experimental Apparatus 

The low-disturbance facility and flared-cone model are first described in this 
section. Then, details of the instrumentation and flow-diagnostic methods employed are 
presented. Lastly, the model mounts designed and built for the present experiments are 
described. 

2.1.1 Low-Disturbance Facility 

The Mach-6 Axisymmetric Quiet Nozzle, which is housed in the Nozzle-Test- 
Chamber Facility (NTC) at NASA Langley Research Center, provided the low-disturbance, 
test environment required for the present stability measurements. A schematic diagram of 
the facility is shown in Fig. 2.1. The locations of the facility fine-screen filter, settling 
chamber, nozzle, and test section are indicated in this figure. More specific details of the 
design and performance of the “quiet” nozzle can be found in Ref. 23. The geometric ratio 
of the nozzle-exit area to the throat area was 56.1. This area ratio yielded a measured Mach 
number of M« = 5.91 in the uniform-flow test region. The stagnation conditions for the 
present tests were Poo© = 130 psia and To«> =810 °R. The Reynolds number per unit length 
for these test conditions was Re/1 = 2.82 x 10 6 ft" 1 . 
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The disturbance environment of the NTC has been documented in a previous 
study. 24 Results of these measurements are relevant to the analysis of the present 
experiments and are summarized in the following. Within the uniform Mach-number 
envelope a large region existed in which flow disturbances were undetectable, as indicated 
by the quiet-flow region shown in Fig. 2.2. Beyond this region, a low-noise test volume 
existed in which disturbances were confined to the 0-50 kHz low-frequency band, 
primarily centered around 16 kHz. The source of these disturbances was acoustic radiation 
from the pre-transitional boundary layer on the nozzle wall. 

2.1.2 Flared-Cone Model 

The NASA-LaRC test model designated the 93-10 flared-cone model was used for 
this study. A schematic diagram of the model geometry is shown in Fig. 2.3. The model 
was constructed with a highly-polished thin skin to minimize effects on the model 
boundary-layer stability due to both surface roughness and heat conduction. The model 
measures 20 H in length and has a 4.6” base diameter. The 10-inch length of the model 
measured from the sharp tip is a straight-walled cone with a 5-degree half angle. The 
remaining 10-inch length of the model is flared outward with a 93.071-inch radius circular 
arc. The flared arc is tangent to the straight wall at the cone-flare juncture. The purpose of 
the flare is to generate an adverse-pressure gradient which may induce transition on the 
model within the operating capabilities of the NTC Facility. 

2.1.3 Model-Surface Instrumentation 

To monitor surface temperature the flared-cone model was instrumented with K- 
type thermocouples. The model-surface thickness directly over the surface-mounted 
thermocouples was 0.03” and 0.06" elsewhere. A total of 51 thermocouples were axially 
spaced along a single meridian at 1-inch intervals between 2 and 9 inches and at 0.25-inch 
intervals from 9 to 19.75 inches, measured from the tip. Along the opposing meridian, 29 
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eq ually -spaced static-pressure taps are located. An electronic data-acquisition unit 
controlled over a standard computer interface was used to monitor and record the 
temperatures. Surface pressures were monitored at select locations on the model with the 
measurement system described in the following. 

2.1.4 Pressure-Measurement System 

To measure the static pressure on the model surface a set of ten MKS 690 Absolute- 
Pressure Transducers was used. The full-scale range of the transducers was 279 lb/ft 2 . 
Stainless-steel tubing which ran through the model sting and outside the test section 
connected the pressure-tap orifices through a manifold system to the transducers. The 
manifold was evacuated by a low-pressure vacuum pump to check the tubing for leaks and 
to provide a vacuum-reference pressure for offset-voltage calibration. This manifold 
system allowed each transducer to be valve-isolated either to the vacuum-reference pressure 
or to its corresponding orifice on the model surface. 

The transducer output voltage was monitored by MKS 670A High- Accuracy Signal 
Conditioners. These signal conditioners allowed the mean-flow surface pressure to be 
monitored with a high degree of precision. The time constant for all measurements was set 
to 400 msec which effectively attenuated frequencies above 0.4 Hz by at least 3 dB. The 
filtered analog signal was then digitally sampled at a rate of F s = 10 Hz. The mean-flow 
pressure level was determined by taking the average of the 100 most-recent samples. This 
provided a 10-second moving average of the surface pressure which mitigated effects of 
random electronic noise in the transducer output and instantaneous pressure fluctuations in 
the physical flow. 

The surface-pressure measurements were accurate to within ± 0.083 lb/ft 2 . The 
freestream static pressure was P« = 13.01 lb/ft 2 , and the freestream dynamic pressure was 
q« = 318.1 lb/ft 2 . From inviscid theory 25 the surface pressure over a straight-walled cone 
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with a 5-degree half angle at a zero-degree angle of attack under these freestream conditions 
is approximately Pcone = 19.74 lb/ft^. Using these values as a reference, the relative 
precision for the non-dimensional pressure coefficient under typical conditions was 
Cp = (Pcone-P«)/q~ = 0-0211 ± 0.00026. 

2.1.5 Schlieren Flow Visualization 

The NTC facility was equipped with a conventional schlieren system for flow 
visualization. The schlieren method visualizes the density variations in a high-speed 
compressible flow. A high-intensity strobe light source was synchronized with a video 
camera which enabled images of the flow to be monitored in real time and recorded on a 
video cassette. The standard frame-capture rate of the video system was 30 Hz. However, 
for each frame the light source was pulsed for only a 1-jisec duration to obtain sharper 
images of the high-speed flow. The schlieren knife edge was oriented parallel to the nozzle 
axis to vis ualize density gradients in the vertical direction. This orientation was chosen to 
highlight the gradients through the model boundary layer in the schlieren image. Because 
most of the model was inside the nozzle, however, the schlieren system was used to 
qualitatively assess the state of the boundary layer over only the aft region. 

2.1.6 Hot-Wire Anemometer System 

The hot-wire anemometer system for obtaining mean-flow and stability 
measurements in the boundary layer consisted of a hot-wire probe, a traverse mechanism, 
and an anemometer circuit. Each of these integral components is described below. 

A schematic diagram of the single-component hot-wire probe used in the present 
study is shown in Fig. 2.4. The tip of the probe consisted of a set of two stainless-steel 
broaches to which the 0.0001-inch diameter platinum-plated tungsten wire was manually 
spot welded. A third broach, as shown in the figure, served as a contact switch which 
completed an electrical circuit when in physical contact with the model. The probe-tip 
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dimensions, dp, dq, and dl in the figure, were manually measured with a precision of 
± 0.001 inch. Dimensions of the probe tips used for the present study are listed in Table 
2.1. The typical length-to-diameter ratio of the thin-filament hot-wire sensors was 200. 

The traverse system positioned the hot-wire probe within the boundary layer on the 
model with a high degree of precision. The interchangeable probe tips were mounted on a 
support arm which was in turn mounted to a set of motor-driven traversing slides. The 
traversing slides were perpendicularly mounted on the ceiling of the test-section box. A 
photograph of the traverse units and the support arm mounted in the NTC test section with 
the 93-10 flared-cone model is shown in Fig. 2.5. This mechanism allowed movement of 
the hot-wire probe tip in a horizontal plane in the test section with a precision of about 
± 0.0004". The traverse system was controlled via a standard computer interface. 

The anemometer circuit used to obtain measurements from the hot-wire sensor was 
a new constant- voltage anemometer (CVA). 2 ** The CVA system was a prototype system 
and remains under continued evaluation and development 27 The selection of the constant- 
voltage system was primarily dictated by its ability to measure disturbances in the 
freestream of the NTC Facility. Other anemometry techniques such as constant-current and 
constant-temperature systems were not capable of detecting the small-amplitude 
disturbances in the freestream of the quiet facility .28 The high-sensitivity prototype CVA 
system was reported to have a 350-kHz bandwidth. 

The CVA output was not calibrated for the present measurements. Interpretation of 
the uncalibrated mean-flow and stability measurements was based on the following 
considerations. Mass flux in the direction normal to the wire and stagnation temperature at 
die sensor both contribute to the anemometer-output signal. The sensitivity of the CVA to 
the mass-flux contribution becomes considerably larger than that to the stagnation 
temperature contribution, however, as electrical heating of the hot-wire sensor increases. 2 ** 
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Electrical heating of the wire is controlled by adjusting the sensor voltage or current The 
present measurements were obtained with the hot-wire sensor operated under high-voltage 
conditions. In addition, the stagnation-temperature fluctuations are assumed to be 
negligible. Thus, the dependence of the anemometer output is primarily attributed to the 
mass-flux contribution. Furthermore, stability theory prescribes a common exponential- 
wave amplification factor for all physical flow-disturbance quantities. Thus, uncalibrated 
measurements of the non-dimensional amplification rates and integrated-growth factors 
may be compared with stability theory independent of the physical flow-disturbance 
quantity. 

Since a single-component hot-wire sensor was used, only the component of flow 
disturbances propagating in the direction normal to the hot-wire itself can be measured. 
Hot-wire measurements were thus restricted to regions where the mean-flow direction was 
known a priori. Therefore, all hot-wire measurements were conducted along the windward 
and leeward meridians, i.e., in the symmetry plane of the axisymmetric model at angle of 
attack. 

2.1.7 Angle-of-Attack Mount 

Since the test facility was not equipped with a mechanism for adjusting the angle of 
attack of the model, a new model mount was fabricated for the present study. The design 
of the new model mount was subject to several constraints which are as follows. The 
traversing mechanism allowed movement of the hot-wire probe only in a horizontal plane, 
parallel to the ceiling of the test section. Thus, the symmetry plane, in which hot-wire 
measurements were to be obtained, was restricted to a horizontal plane. In addition, the 
hot-wire probe-arm mount and the software algorithm which controlled the traverse 
movement were configured to conduct measurements only on the left side of the model 
when looking upstream into the nozzle. So, to obtain measurements on both windward 
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and leeward rays, the axisymmetric model was physically required to be placed at both 
positive an d negative side-slip angles in the test section to investigate the effects of angle of 
attack. 

Schlieren photography of the windward and leeward meridians could not be 
obtained concurrent with hot-wire measurements when the symmetry plane was horizontal. 
In this configuration, the light path of the NTC schlieren system through the horizontal 
symmetry plane was obstructed by the model. Therefore, in order to obtain schlieren 
visualization the model was also required to be placed with the windward and leeward rays 
in a vertical plane. As an additional constraint, the model was required to remain within the 
uniform Mach-number envelope of the nozzle and upstream of the nozzle-exit shock. 

The maximum angle of attack attainable, subject to the above constraints, was 
determined to be 5 degrees. The pivot point was located on the nozzle axis 4-Vi6" 
upstream of the nozzle exit. Thus, the following angle-of- attack test cases were 
established: a = 0°, 2°, and 4°. The new model mount provided for the positive and 
negative physical sideslip angle configurations, {3 = i 2° and i 4°, for hot-wire 
measurements. Additionally, mounts were needed for the physical angle- of- attack 
configurations, a = 2° and 4°, for schlieren visualization. Here, a positive angle denotes 
the windward case, and the negative angle refers to the leeward case for hot-wire 
measurements. An existing mount was used for the zero-degree baseline case. 

The above requirements were satisfied by the design of a set of mounting blocks. 

A common base block was mounted to the floor of the test section, and interchangeable 
blocks were bolted to the base block for sideslip and angle- of- attack configurations. An 
existing strut mounted the model to the new blocks. For hot-wire measurements, one top 
block allowed the model to be mounted at sideslip angles of [3 = +2° and -4° as well as |3 =- 
2° and +4° when inverted. Two other top blocks were machined with mounting surfaces at 


14 



angles of attack of a = 2° and 4° to be used for schlieren visualization. A photograph of the 
short strut and the four new model mounting blocks is shown in Fig. 2.6. The blocks 
were precision machined from high-strength Aluminum alloy, and steel-threaded inserts 
were used in all bolt holes. Three-view engineering drawings of the blocks are included in 
section 8. 1 of the Appendices. The typical precision of the model angle with the new 
mounting blocks was about ± 0.2°, equivalent to the precision of the existing zero-degree 
mount A photograph of the flared-cone model mounted on the new blocks in the NTC 
Facility is shown in Fig. 2.7. 

2.2 Experimental Procedure 

The model was tested at angles of attack of a = 0°, ± 2 °, and ±4°. Baseline 
measurements were obtained along a single meridian for the a = 0° case. Angle-of-attack 
measurements were obtained along the windward and leeward meridians for the cl — 2 and 
4° cases. Surface and boundary-layer measurements characterized the mean-flow state. 
Also, boundary-layer stability measurements documented the development of flow 
disturbances. Table 2.2 lists the specific measurements obtained with the above 
diagn ostics for each test case. Procedures followed to conduct the experiments are outlined 

below. 

2.2.1 Mean-Flow Surface Measurements 

Mean-flow surface-temperature measurements verified the thermal-equilibrium state 
of the model. Also, mean-flow surface-pressure measurements documented the adverse- 

pressure gradient on the flared-cone model. 

To determine the thermal-equilibrium state, the surface temperature was monitored 
during a wind-tunnel run at test conditions until the temperatures remained constant with 
time. Prior to starting hypersonic flow in the tunnel, a low-speed preheat cycle with air- 
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supply conditions of Po = 50 psia and To = 810 °R, which was routinely applied to alleviate 
condensation, was then used to bring the model to the documented equilibrium 
temperature. Preheating the facility insured that the model surface would quickly reach a 
thermal-equilibrium state, thus maximizing data-acquisition time during hypersonic-flow 
operation. All measurements presented herein were performed with the model in a thermal- 
equilibrium state. 

Prior to conducting pressure measurements, each signal conditioner was re- 
calibrated to compensate for voltage drift and variations in ambient conditions. The 
transducers were initially valve-isolated to the evacuated manifold during the preheat cycle 
and the start of hypersonic flow. After stable hypersonic flow was established, the isolator 
valves were then turned to establish the hydraulic connection from each transducer to the 
corresponding pressure orifice. Time-averaged surface-pressure measurements were 
monitored until these reached steady-state values. 

2.2.2 Mean-Flow Boundary-Layer Measurements 

Off-surface measurements were obtained with the hot-wire anemometer system. 
Time-averaged mean and rms signals were recorded from the CVA system. The mean 
output qu alitativ ely represented the mean-flow mass-flux profile. The rms output is a 
measure of the flow-disturbance energy. Thus, the peak-rms location identified the locus 
of maximum disturbance energy. To optimize the SNR of the disturbances, subsequent 
stability measurements were obtained at these maximum-energy locations. 

Precision movement of the hot-wire probe and data acquisition from the CVA 
system were controlled over a standard computer interface. For the baseline and leeward 
cases the probe was traversed over 17 stream wise locations in 0.5” increments along the 
nozzle axis. Fewer stations were surveyed for the windward cases because of the confined 
space between the model and the nozzle wall. Also, the location of the hot-wire 


16 



measurements was restricted to the flared region of the model for all test cases. At each 
streamwise location the CVA output was recorded at 13 discrete locations within the 
boundary layer. The distribution of these discrete points was smoothly clustered near the 
peak-rms location. The peak-rms location was first determined by a short survey over the 
most-downstream stations. A time delay was applied after any incremental probe 
movement and before any hot-wire measurements were recorded. This delay insured that 
transient effects of electromagnetic interference from the traverse drive motors did not 
corrupt the measurements. 

2.2.3 Boundary-Layer Stability Measurements 

The hot-wire anemometer system was also used to measure stability characteristics 
of the boundary layer. These stability measurements were based on time histories of the 
instantaneous flow-disturbance fluctuations at the maximum-energy locations. The time- 
series data were recorded as follows. The CVA output was first analog filtered with a 
pass-band from 1 kHz to 810 kHz to minimize the effects of low-frequency electronic noise 
and high-frequency aliasing. Then, 40,000 data values were digitally sampled at a rate of 
F s = 2 MHz with 8-bit numerical precision with an oscilloscope. 

2.3 Navier-Stokes Computations 

The laminar flow over the flared-cone model at angle of attack was modeled by the 
numerical solution of the Navier-Stokes equations. These computations were performed 
on the North Carolina Supercomputing Center’s Cray Y-MP and Cray T90. The following 
sections present details of the flow solver, the grid and boundary conditions, and the test 

cases computed. 
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2.3.1 Flow Solver 


The NASA-developed code CFL3D was used in the present study. 29-31 CFL3D 
has been used to obtain mean-flow solutions on which LST computations were based for 
the flaied-cone model at a zero-degree angle of attack in an earlier study. 32 Also, the 
suitability of CFL3D for modeling general high-speed flows has been previously 
assessed. 33 Salient features of the code are now discussed. 

Steady-state solutions were obtained with local time stepping. A flux-difference 
splitting scheme was applied for spatial differencing of the Euler fluxes with upwind-biased 
third-order differencing and smooth flux limiting. Laminar viscous fluxes were included in 
the three coordinate directions simultaneously. The scalar matrix-inversion method was 
applied. Convergence-acceleration techniques such as multi-gridding and mesh sequencing 
were also utilized. 

2.3.2 Grid & Boundary Conditions 

An algebraic mesh generator was used to compute coordinates for the half-body 
plane-symmetric grid. A hidden-line plot of the grid is shown in Fig. 2.8. The farfield 
boundary is shown reflected about the X-Z plane for clarity. Grid points were smoothly 
clustered near the tip of the model and near the surface. Circumferential angular spacing of 
planes common to the centerline was held constant Freestream conditions were prescribed 
on the inflow and farfield boundaries. Plane-symmetry boundary conditions were applied 
to the leeward and windward planes. Over the model surface, adiabatic and no-slip 
boundary conditions were imposed. At the outflow plane a one-point extrapolation was 
specified. 

Previous computations have demonstrated the application of CFL3D to model the 
mean-flow over a 5-degree half-angle cone at a 2-degree angle of attack under similar 
freestream conditions.^ 4 Results of the grid-refinement study conducted for this previous 
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investigation were used as a guide in choosing grid dimensions for the present study. The 
present grid dimensio ns were 109 x 41 x 105 in the stream wise, circumferential, and 
surface-normal directions, respectively. These dimensions allowed for two coarser meshes 
to be used in the mesh-sequencing process and in multi-grid acceleration at the finest mesh 
level. 

2.3.3 Test Cases & Freestream Conditions 

The flow was computed for a = 0°, 2°, and 4° for comparison with corresponding 
experiments. Additionally, to assess the impact of the limited precision of the model 
alignment, the case where a = 0.2° was also computed. The prescribed freestream 
conditions were M« = 5.91, To* = 101.9° R, and Re/1 = 2.82 x 106 ftr* which 
corresponded to the wind-tunnel test conditions. The solutions were advanced until the 
adiabatic surface-temperature distributions converged to a specified tolerance. 

2.4 Experimental & Computational Data Processing 

This section describes the post processing applied to the experimental and 
computational (fa ta, Transformation of the hot-wire coordinates and evaluation of the 
boundary-layer thickness and the second-mode disturbance frequency are discussed in the 
following sections. Estimation of flow-disturbance power spectra and bispectra is then 
discussed. Finally, this chapter is concluded with descriptions of the computation of flow- 
disturbance amplification rates and integrated-growth factors. 

2.4.1 Hot-Wire Coordinate Transformation 

The contact switch was used as a coordinate reference for the hot-wire 
measurements. However, the hot-wire measurement location was not coincident with the 
contact-switch location. Also, these coordinates were measured with respect to the nozzle 
coordinate system. Thus, a transformation was necessary to reference the measurements to 
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model arc-length and surface-normal coordinates. The algorithm of this coordinate 
transformation is described in section 8.2 of the Appendices. 

The contact-switch coordinates were reliable within the precision of the traverse 
mechanism, about ± 0.0004 inch. However, the coordinate transformation required the 
use of limited-precision, manually-measured dimensions, and these are as follows. The 
axial position of the model relative to the nozzle-exit plane, (Xbase~Xexit)» was measured 
with a precision of ± 0.0156 inch, and the sideslip angle, p, was accurate to within ± 0.2°. 
The angle of the hot-wire, probe arm, 0, was adjusted with an accuracy of ± 0.5°. And, 
dimensions of the hot-wire probe tips, dp and dq, were measured to within ± 0.001 inch. 

The accuracy of the transformation was estimated by adjusting the values of the 
manually-measured dimensions over the range of their respective tolerances and examining 
the net effect on the final coordinates. The reliability of the transformed surface-normal 
coordinates was estimated to be about i 0.002 inch and the accuracy of the surface, arc- 
length distances was estimated at ± 0.02 inch. As a reference scale, the boundary-layer 
thickness near the base of the 20-inch long model at a zero-degree angle of attack is about 
0.06 inch. 

2.4.2 Boundary-Layer Thickness Estimation 

The boundary-layer thickness, 5, was independently estimated from both mean- 
flow measurements and Navier- Stokes computations. The measured estimate was 
determined by the surface-normal distance where the CVA, mean-output profile approached 
a near-zero slope, to within a specified tolerance. The mean output was primarily attributed 
to mass flux. Because the discrete measurement locations were clustered near the peak-rms 
location, the data near the boundary-layer edge were sparsely spaced. Application of a 
smoothing cubic-spline fit determined 5 with higher resolution than the spacing of the 
original data. However, the smoothing fit tended to bias the edge location outward by a 
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small amount Thus, the reliability of the measured 8 estimate was limited by the accuracy 
of both the coordinate transformation and the curve fit the combined effects of which 
resulted in a net accuracy of about ± 0.005 inch. 

The computational estimate was evaluated from the computed boundary-layer 
velocity profile. Evaluating 8 based on the slope of the velocity profile was a practical 
method because the result was not highly sensitive to the edge-gradient tolerance. While a 


velocity gradient existed between the boundary-layer edge and the model shock, this 
gradient is small when compared to that within the shear layer itself. The variation in 8 due 
to adjusting the edge-gradient tolerance fell well within the precision of the experimental 


measurements. An interpolating cubic-spline fit was applied to the numerical data to further 
increase the resolution of the 8 estimate beyond that of the computational mesh. 

2.4.3 Second-Mode Disturbance-Frequency Estimation 

For a general wave-like disturbance, the frequency, F, and the wavelength, X, are 


related to the speed of propagation, IV1, by 

IVl = X F. (2-1) 

The second-mode disturbance wavelength has been shown to scale with the 


boundary-layer thickness according to *6, 17, 32 

X2ndMode ss 28. (2.2) 

The speed of propagation of the second-mode disturbance was estimated by the 
speed at the location of maximum flow-disturbance energy within the boundary layer. The 
boundary-layer disturbances are assumed to be at a maximum amplitude and energy at the 
most-unstable location in the boundary layer. This location is referred to as the 
generalized-inflection pointA 35 The generalized-inflection point, ygip, is the location 
where 
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(2.3) 


a / aivA 
dy\ ay) 


= o. 


Combining Eqs. 2.1 and 2.2 above, the frequency of a second-mode disturbance, 
F?nfi Mode» can be estimated based on the velocity magnitude at ygip, IV gipl, and the 


boundary-layer thickness, 8, by, 

* Wgip 1 n 

Fmd Mode * — “ . ( 2 - 4 ) 

28 

Row solutions obtained from CFL3D along the windward and leeward meridians 
were used to evaluate IVgipI and 28. The above estimate is not presented as a substitute for 
more-detailed stability computations. Without LST-based computations available for 
comparison with the present angle-of-attack test cases, however, this estimate provided a 
general indication of the range of frequencies that may be measured in the experiments. 
Moreover, this result provided a basis for the assessment of angle-of-attack effects on the 
second-mode disturbance frequency. 

2.4.4 Power-Spectrum Estimation 

To identify flow-disturbance modes, averaged-power spectra of the time-series, 
stability measurements were computed as follows. The 40,000 sample record at each 
measurement location was divided into segments of 256 samples. This segment length 
corresponds to a duration of T = 0.128 msec for the 2-MHz sampling rate. From each 
segment the mean value of the segment was subtracted. A Hanning window was then 
applied to the segment The fast-Fourier transform of the windowed segment was 
computed with 1024 data points. The resulting complex-valued spectrum was multiplied 
by its conjugate to yield the real-valued power spectrum. Power spectra from 156 non- 
overlapping segments of the original record were averaged to compute the averaged power 
spectrum. This procedure was applied to an additional record of 40,000 points recorded at 
the same location, and these spectra were then averaged for each measurement location. 
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The frequency-sampling interval was AF = 1.95 kHz for the computed spectra. 

The frequency resolution of the Hanning window was 4 H ~ 3 1.25 kHz, or twice that of a 
rectangular window of the same length.^ However, die Hanning window provided more 
than twice the amount of side-lobe attenuation. Thus, the effective dynamic range was 
increased with the use of a Hanning window over the rectangular window. Such a feature 
was useful in detecting small-amplitude disturbances which may have been otherwise 
obscured by frequency leakage from stronger signal components. 

2.4.5 Bispectrum Estimation 

The bispectrum was used to detect non-linear phase-coupling disturbance 
interactions. Phase coupling occurs when the sum of the phases of frequency components 
Fi and F 2 is equivalent to the phase of the sum frequency, (Fi + F 2 ). This phase 
relationship is described by 

Z X(F X ) + Z X(F 2 ) = Z X(Fi+F 2 ), (2.5) 

where X(F) is the Fourier transform of the time-series signal x(n). Unlike the first-order 
power spectrum, the higher-order bispectrum retains phase information about frequency 
components as shown below. The bispectrum of x(n) is defined as 

B(Fi,F 2 ) = E[X(Fi) X(F 2 ) X*(Fi+F 2 )], (2.6) 

where E[ ] indicates the expected-value statistical operation. Bispectra were computed 
using an fft-based algorithm similar to the power-spectrum estimation method. The time- 
series data record was divided into segments, and for each segment X(F) was evaluated. 
The triple product, [X(Fi) X(F 2 ) X*(Fi+F 2 )L was computed for each segment and then 
averaged. The averaged estimate for B(Fi,F 2 ) was obtained from 19 segments of 256 
points each. The frequency-sampling interval was AF = 1.95 kHz and the frequency 
resolution was 16.125 kHz. However, the bispectrum estimate was smoothed in the 
frequency domain with a square convolution mask approximately 8 kHz wide, so the net 
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resolution was about 25 kHz. This resolution was sufficient to identify frequency bands of 
phase coupling. 

2.4.6 Amplification Rates 

The boundary-layer stability was also characterized by the local flow-disturbance 


amplification rate. 


In its non-dimensional form the amplification rate, -otj, is given by 


-<Xi = 


jl3A 


(2.7) 


where A is the amplitude of the flow disturbance and R is the stability Reynolds number 
defined by R = (Re s ) 0 - 5 . The above non-dimensional form is independent of any locally- 
constant multiplicative calibration factor, however, additive noise may obscure detection of 
flow-disturbance amplification. 

The amplitude spectra were first evaluated by taking the square root of the average- 
power spectra. A smoothing cubic-spline fit was then applied to refine the data along the 
streamwise direction from 17 to 64 locations for each discrete frequency. A second-order 
accurate central finite difference was used to approximate the spatial derivative, dAIdR, 
from the fitted amplitude data. The central values of the amplitude were used for the 
normalization coefficient, (V2A)* Using j to indicate the array index along the streamwise 


direction the above finite-difference approximation can be written as 

( x 1 ( A j+1 - A J-l) 

( l)j 2Aj (Rj+i - Rj-i) 


2.4.7 Integrated-Growth Factors 


( 2 . 8 ) 


The flow-disturbance integrated-growth factor, commonly referred to as the N 
factor, locally measures the net growth of a flow disturbance from the initial point of 
amplification. Detection of the initial growth of infinitesimal disturbances is limited by the 
finite level of noise in the CVA system. Thus, the absolute growth factor is not physically 


measurable. Using the most-upstream location as a reference, however, the relative growth 
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factor at downstream locations can be measured. The N-factor difference was evaluated 
from the measured amplitude spectra according to the relationship shown previously in Eq. 
1.3 and repeated here with the following change in notation. N is the local integrated- 
growth factor, and No is the unknown integrated-growth factor at the arbitrary reference 

location, xo. Their difference is given by 

N - N 0 = In (2.9) 

Ao 

where A is the local amplitude, and Ao is the amplitude at the reference location, typically 
chosen as the most-upstream measurement location. To experimentally detect growth the 
disturbance must have an amplitude greater than the additive noise in the anemometer 
system. Thus, regions where growth is not experimentally measured should not 
necessarily be interpreted as regions of neutral stability. Downstream of measured 
disturbance growth, the slope of the N-factor difference can be compared on a consistent 
basis with LST predictions. 

In regions where the flow signal is experimentally detectable above the noise floor, 
the measured N-factor difference is not degraded by the spatially-invariant additive noise. 
Electronic noise is assumed to be independent of streamwise location, so (N-No)noise 
would be negligible. However, the flow-disturbance environment of the test facility is not 
spatially invariant Consequently, the impact of facility noise on the present stability 
measurements was assessed by monitoring disturbances in the 0-50 kHz low-frequency 
band. 
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3 Results & Discussion 


For all angle-of-attack test cases the freestream unit-Reynolds number was Re/1 = 
2.82 x 10 6 fr 1 . Under these conditions the spatial distribution of the stability Reynolds 
number, R, is plotted versus the arc-length Reynolds number, Re s , with the flared-cone 
model geometry in Fig. 3.1. The cone-flare tangent point at X = 10" corresponds to Re$ — 
2.35 x 10 6 , or R = 1533. Surface measurements were obtained over the length of the 
model. Hot-wire measurements, however, were conducted over the flared region typically 
in the range X =1 1"-19" which corresponds to (Re s /10 6 ) = 2.6-4.5 or R =1600-2100. All 
experimental and computational data presented herein are plotted against (Res / 10^) to 
indicate streamwise location. 

The discussion of results below is organized as follows. First, Navier-Stokes 
computations of the flow field over the flared-cone model at angle of attack are presented. 
Next, surface measurements and computations which document the mean-flow state of the 
boundary layer are presented. Schlieren flow-visualization images are then presented 
which verify the boundary-layer state. Following this, boundary-layer measurements 
which characterize the mean flow and the stability are discussed. Finally, the effects of 
angle of attack on the boundary-layer stability and transition for the flared-cone model are 
summarized. 

3.1 Flow-Field Computations 

Navier-Stokes computations are presented here to demonstrate the unique features 
of the flow over the flared-cone model. Also, the general effects of angle of attack on the 
mean flow are discussed. A flooded-contour plot of the computed pressure field is shown 
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for the zero-degree case in Fig. 3.2. The color map along the top of the figure indicates the 
numerical values of the non-dimensional pressure coefficient Contour lines of constant 
pressure are shown on the surface. Axial planes are shown at a station 5 from the tip and 
at the outflow plane. The solutions in the windward and leeward planes are identical for 
this axisymmetric case. The upstream portion of the model is equivalent to a straight- 
walled cone and exhibits similar features. The pressure rises across the shock originating 
from the tip. Pressure contours follow rays emanating from the tip, including the cone 
surface, consistent with inviscid-cone theory.^ Downstream of the tangent, however, the 
pressure rises at a discontinuous rate due to the adverse-pressure gradient induced by the 
flare's outward curvature. Compression waves propagate away from the flare surface. A 
slight degree of asymmetry is observed for the a = 0.2° case as shown in Fig. 3.3 . The 
skewed surface-pressure contours indicate that the pressure increases on the windward ray 
and decreases on the leeward ray relative to the axisymmetric zero-degree case. 

Computations for the a = 2° case are shown in Fig. 3.4. The pressure has 
increased substantially on the windward ray and decreased on the leeward ray compared to 
the zero-degree case. Pressure contours in the leeward plane are perpendicular to the 
surface because the pressure gradient in the surface-normal direction is very small within 
the viscous boundary layer. The vertical extent of these normal pressure contours indicate 
the approximate thickness of the boundary layer over the flare. This effect was not clearly 
shown in the baseline computations because the boundary layer was much thinner. The 
boundary layer on the leeward ray grows even thicker in the a = 4° case as shown in Fig. 
3.5. Conversely, the increased pressure levels on the windward ray result in a much 
thinner boundary layer although this feature is not clearly shown in these plots. The 
discrete mesh fully contains the conical shock for all of the angle- of- attack cases 
investigated. The pressure field on the flared-cone at angle of attack causes the flow to 
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migrate from regions of high pressure towards regions of low pressure. The resulting flow 
asymmetry may significantly influence the development of instabilities in the boundary 
layer. 

3.2 Mean-Flow Surface Measurements & Computations 

The mean-flow surface measurements were obtained with the use of the model's 
internal instrumentation. Surface pressure and temperature distributions for the flared-cone 
at angle of attack are presented in the following sections. 

3.2.1 Pressure Distributions 

The measured surface-pressure distribution for the oc = 0° baseline case is plotted in 
Fig. 3.6. Overlaid on this plot are the Navier-Stokes computations for the a = 0° and 0.2° 
cases. Pressures were m easured at two locations on the constant-pressure cone region. 

The upstream data point falls close to the a = 0.2° leeward computations; conversely, the 
other data point falls above the a = 0.2° windward computations. This contraindication 
illustrates the difficulty encountered in detennining actual model alignment from surface- 
pressure data. Because each pressure is measured independently, however, unique 
c ha racteristics of the hydraulic connections and transducers can result in a significant bias, 
as observed in the present data. Over the flare the measured distribution shows a slightly 
higher pressure rise than the laminar computations. This discrepancy may be due in part to 

transition onset or the physical alignment of the model. 

The oc = 2° pressure distributions are plotted in Fig. 3.7 . The measurements show 
good agreement with the computations along the cone portion. The measured pressure 
along the windward ray falls slightly higher than the computed values. The measured 
leeward dis tribution shows a more significant deviation from the computed laminar trends, 
especially over the flare. Towards the base these measurements indicate a larger adverse- 
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pressure gradient than that predicted by the computations. This larger adverse-pressure 
gradient is due to transition onset Because of the sparseness of the data the location of 
transition onset could not be verified with high precision. Surface-pressure measurements 
were not obtained for the a = 4° configuration. 

3.2.2 Temperature Distributions 

The surface-temperature distributions for cc = 0°, cc = 2°, and cc = 4° are plotted in 
Figs. 3.8, 3.9, and 3.10, respectively. The surface temperature, T, is normalized by the 
freestream stagnation temperature, Tooo. The vertical interval between the (TATooo) grid 
lines corresponds to a temperature change of AT = 16° R. Gradients in the computed 
distributions near the tip are artifacts of the singularity in the numerical mesh and are not 
physically significant 32 The ± 2° R precision of the temperature measurements 
corresponds to ± 0.0025 x To~ for the normalized temperature, which is equivalent to the 
height of the plot symbols. 

Navier-Stokes computations for the a = 0° baseline case, shown in Fig. 3.8, 
exhibit the following features. Over the cone portion of the model the temperature is fairly 
constant At the cone-flare tangent the slope of the temperature distribution increases 
discontinuously. Over the flare, the temperature increases moderately by about 1.4%. The 
estimated variation due to the model alignment is indicated by the envelope of the <x = 0.2 
computations. The largest variance occurs at the tangent where the windward distribution 
falls 0.5% higher, and the leeward distribution falls 0.4% lower than the baseline case. 
The a = 0° measurements fall within this envelope up to Res = 3.8 x 1(A Downstream of 
this location transition onset is indicated by the large temperature increase. The transition- 
onset location is determined to be Res = 4.2 x 10^ by the intersection of the asymptotes 
extended from the laminar and transitional regions. The development of fully-turbulent 
flow represented by a smooth asymptotic distribution downstream of the transitional peak 
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is not observed, however. The temperature decrease at the last measurement station may 
indicate that transition is nearing completion. Near the base, however, conductive cooling 
from the internal model structure may also produce a similar effect 

The a = 2° distributions are plotted in Fig. 3.9. Over the cone portion the 
computations predict a greater temperature difference between the windward and leeward 
rays than that observed in the experiment In this region, however, the physical 
dimensions of the model are small and the windward and leeward rays are strongly 
conductively coupled. Such thermal conduction tends to average out temperature 
differences between the windward and the leeward rays near the tip. The computed 
windward temperature is essentially constant over the cone and increases only 0.8% over 
the flare. The windward measurements fall well below the computations near the tip, but 
asymptotically approach the computations near the base as thermal conduction with the 
leeward ray is mitigated. The small rise in the measurements near the base may be due 
either to transition onset or to conduction effects. In contrast, the leeward computations 
show a sharp rise in the surface temperature at the tangent suggesting laminar separation. 
Separation is not observed in the experiment, however. Conduction effects may somewhat 
smooth out the distribution, but not to the extent that such a separated region would be 
undetectable. A more likely explanation is that the boundary layer remains attached as a 
result of the high-energy flow associated with the pre-transitional boundary layer. 
Transition onset is detected from the measured distribution at Re s = 3.33 x 10 6 . 

Fig. 3.10 shows the a = 4° results. Strong conduction effects are again observed 
as the windward and leeward temperatures have equilibrated near the tip. The windward 
measurements follow laminar trends but are shifted by conduction effects. A slight 
temperature rise near the base may indicate either transition onset or conduction. The 
leeward measurements are relatively constant along the cone and sharply rise at the tangent 
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Just downstream of the tangent, at Reg = 2.75 x 10 6 , separation is indicated by the region 
of constant temperature. The leeward computations predict separation slightly upstream of 
the measured separated region. Separation acts a bypass mechanism resulting in 
instantaneous transition to turbulence for this test case. 

Closer examination of transition onset for the cc = 0° baseline and ot = 2° leeward 
cases reveals a smooth increase where transition onset occurs. In contrast, the ot = 2 and 
4° windward cases both exhibit a more discontinuous temperature rise near the base. This 
observation supports the conclusion that the windward cases remain fully laminar, and 
conduction effects are responsible for these observed temperature rises. In addition, the 
laminar computations predict separated regions for both leeward cases. However, 
separation was detected only for the ct = 4° leeward case. Conduction effects were not 
strong enough to completely dissipate the constant-temperature region in the distribution for 
this case. Thus, such a region would also be detectable, if one indeed existed, for the ex — 
2° leeward case. Based on these considerations, one may conclude that the transitional 
boundary layer remains attached for the a = 2° leeward test case. 

3.3 Schlieren Flow Visualization 

A schlieren image of the flow over the base region of the model is shown in Fig. 
3.11 for the a = 0° baseline case. In this schlieren image the flow is from the left to the 
right The nozzle exit is on the left of the image and the model surface is along the bottom 
of the image. The boundary layer appears as a dark band very close to the model surface. 
Additional intensity variations are due to the turbulent shear layer and nozzle-exit shock 
surro undin g the model as well as imperfections in the test-section windows. These 
background intensity variations are distinct from the boundary layer on the model. Near 
the nozzle-exit plane the edge of the boundary layer is clearly smooth. The local Reynolds 
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number in the nozzle-exit plane was Re s = 3.8 x 10 6 which was within the laminar region. 
The observed transition location of Res = 4.2 x 10 6 corresponds roughly to the middle of 
the field of view. Intermittent wave-like structures are visible in the transitional boundary 
layer closer to the base. The wavelength of these structures is approximately twice the 
boundary-layer thickness; thus, these boundary-layer disturbances are attributed to the 
second-mode instability. The end of transition and the establishment of fully-turbulent 
flow is not apparent in this schlieren image. 

Transition was not observed on the windward ray at angle of attack. The schlieren 
images shown in Figs. 3.12 and 3.13 verify die laminar state of the boundary layer on the 
windward rays for the a = 2° and 4° test cases, respectively. The field of view of the 
schlieren system was adjusted to visualize die flow over the windward ray located on the 
bottom of the modeL In these images the nozzle exit is located on the left, and the model is 
now along the top of the image. The boundary layer appears as a bright smooth band along 
the bottom of the model The nozzle-exit shock and shear layer again contribute to 
background intensity variations. Also, the shock induced by the model can be faintiy seen 
in the lower left as a dark line almost parallel to the model surface. This verifies that the 
shock did not reflect inside the test nozzle and impinge on the model surface. The nozzle- 
exit shock approaches the base of the model as the angle of attack increases. At a 4-degree 
angle of attack, the shock just misses the base of the model. 

Transition onset was observed well upstream of the nozzle-exit on the leeward ray 
for the a = 2° case. The schlieren image for this turbulent case is shown in Fig. 3.14. In 
this image the boundary layer on the model is not clearly discernible from the background 
intensity variations. Closer inspection of the image near the cone surface reveals 
intermittent structures that are slightly more prominent than the background intensity 
variations. These features may be attributed to turbulent structures in the boundary layer. 
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The turbulent boundary layer on the leeward ray for the 4-degree case is also difficult to 
distinguish from background intensity variations in the schlieren image. Thus, the image 
for this test case is not shown here. 

3.4 Mean-Flow Boundary-Layer Measurements & 
Computations 

Time-averaged measurements obtained within the boundary layer consist of both 
the mean and rms output from the CVA system. The mean profiles provided a basis for 
estimating the boundary-layer thickness. Also, the rms output is a measure of the flow- 
disturbance energy. The uncalibrated C VA-output voltages were arbitrarily normalized and 
scaled for the mean and rms profile plots presented below. 

3.4.1 Boundary-Layer Thickness 

The mean output is plotted in Fig. 3.15 against the surface-normal distance, 
indicated by Y, at several streamwise locations for the a = 0° baseline case. The profile 
shape, which is primarily dependent on mass flux, does not appreciably change in the 
t fm-iinar region. Downstream of transition onset, Re$ = 4.2 x 10^ for this case, the mean- 
profile shape broadens. The measured 5 distribution is shown in the Y-Re s plane of this 
3-D plot. The 5 measurements are plotted in Fig. 3.16 with Navier-Stokes computations 
for the a = 0° and 0.2° cases. The computations show conical boundary-layer growth up to 
the tangent The adverse-pressure gradient retards boundary-layer growth over the flare, 
however. The measurements fall within the envelope of the ct = 0.2° computations over the 
most of the flare. Also, the a = 0° computations fall within the precision of the 
measurements indicated by error bars in the plot The experimental data do fall slightly 
below the computations at the upstream locations. This discrepancy indicates that the 
model was aligned at a slight windward angle. Additionally, transition onset is indicated 
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by the increase in 8 near Re s = 4.2 x 10 6 , contrary to the decreasing trend of the laminar 
computations. 

Mean profiles for the a = 2° and 4° windward cases are shown in Figs. 3.17 and 
3.18, respectively. Both of these cases follow similar trends. The profile shape does not 
appreciably change in the stream wise direction. Also, the measured 8 estimate remains 
fairly constant with the exception of some experimental scatter. The 8 measurements are 
plotted with Navier-Stokes computations in Fig. 3.19. The increased pressure level on the 
windward ray further suppresses boundary-layer growth over the cone relative to the 
baseline case. Over the flare 8 remains constant, as opposed to the decreasing trend in the 
baseline case. These measurements do not show any indication of transition onset for 
either windward case. The measured 8 falls slightly higher than the computations. This 
outward bias is due in part to the sparseness of the raw data and the application of the 
smoothing-spline fit 

The a = 2° and 4° leeward mean profiles are shown in Figs. 3.20 and 3.21 , 
respectively. These cases exhibit the following trends. The upstream profiles have a 
laminar shape with a well-defined boundary-layer edge. However, as the profile evolves 
downstream the shape fills out and the edge falls beyond the Y-range over which hot-wire 
measurements were obtained. The fuller shape represents increased mass-flux levels closer 
to the surface which are characteristic of turbulent boundary layers. 

3.4.2 Disturbance-Energy Profiles 

The disturbance-energy profiles, represented by the CVA-rms output, for the cc = 0° 
baseline case are plotted in Fig. 3.22. These profiles show a localized disturbance 
developing within the boundary layer. Downstream of transition onset, Re$ = 4.2 x 10^ 
for this case, the localized disturbance grows in energy and broadens in extent The Y 
locations where flow disturbances have a peak amplitude, indicated by the circular plot 
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symbols at the rms maxima in the figure, correspond to the most-unstable locations in the 
boundary layer. To illustrate this correspondence the loci of the rms maxima are plotted 
with the generalized-inflection point, Ygip, distributions for the a = 0° and 0.2° 
computations in Fig. 3.23. The computed distributions are somewhat affected by errors 
due to interpolation from the discrete mesh, but remain within the precision of the 
measurements. 

Rms profiles for the a = 2° and 4° windward cases are shown in Figs. 3.24 and 
3.25, respectively. In the a = 2° case, the rms profiles have two maxima. The peak closer 
to the boundary-layer edge is attributed to the disturbances at the generalized-inflection 
point The peak closer to the surface results from the increased sensitivity of the hot-wire 
sensor under decreased mass-flux conditions. This secondary peak is also apparent to a 
lesser degree in the a = 4° windward case. Both cases show that the disturbance energy is 
confined to a small region within the boundary layer. As previously documented by mean 
profiles and surface measurements, flow disturbances have not grown sufficiently to 
significantly alter the laminar mean flow for the windward cases. 

In contrast with the windward laminar cases, the rms profiles for the a = 2° and 4° 
leeward cases are shown in Figs. 3.26 and 3.27, respectively. The a = 2° case illustrates 
the evolution of the flow disturbances through the transition process. Initially laminar, the 
localized disturbance is confined to the vicinity of the generalized-inflection point This 
localized disturbance quickly spreads within the boundary layer, however, and grows 
substantially in the downstream direction. The upstream profile for the a = 4° case shows 
a more widespread disturbance, which is expected for this turbulent case. 
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3.5 Boundary-Layer Stability Measurements 

The stability measurements obtained from the time-series hot-wire data are 
presented below. First, the a = 0° baseline case is presented, followed by the windward 
and leeward test cases. The flow-disturbance power spectra at the maximum-energy 
locations are examined at several stream wise stations to determine the development of 
individual disturbance modes. Also, facility noise in the 0-50 kHz band is monitored to 
determine what effect these freestream disturbances may have on the stability 
measurements. The power spectra of the electronic noise in the CVA system is plotted in 
Fig. 3.28. This shows the typical shape of the noise floor above which flow-disturbances 
must rise to be detectable. 

3.5.1 a ss 0° Baseline Case 

Power spectra for this case are plotted in Fig. 3.29. The normalized power, P/Po, 
is plotted on a decibel scale along the vertical axis against the frequency, F, for each 
stream wise measurement location, Res. At the most-upstream location flow disturbances 
are small, and the measured signal consists primarily of electronic noise. Downstream, 
boundary-layer disturbances are detected around 266 kHz. Another significant disturbance 
falls within the 0-50 kHz band of the facility noise centered around 15.6 kHz. The 
normalized amplitude spectra, ln(A/Ao), are plotted in Fig. 3.30. Plotting the spectra in 
this form shows more clearly growth in the frequency band of the sub harmonic of the 266- 
kHz disturbance, around 133 kHz. Also, disturbance growth is detected for higher-order 
harmonics of the 266 kHz disturbance at 531 kHz and 797 kHz. The appearance of these 
frequency-coupled disturbances suggests the onset of non-linear disturbance interactions. 
Detection of these higher-frequency disturbances is unexpected, however, for the following 
reasons. Not only do these harmonics fall well beyond the reported frequency response of 
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the CVA system, 26 but also these unstable disturbances are not predicted by LST 
computations. 32 

The second-mode frequency estimates based on the present Navier-Stokes 
computations are overlaid on an image of the ln(A/Ao) data in Fig. 3.3 1. The frequency is 
along the vertical axis, and the location, indicated by Re s , is along the horizontal axis. 

Both the a = 0° baseline and the a = 0.2° windward and leeward computations are shown 
with the a = 0° measurements. The levels of measured growth are indicated by the color 
map on the right of the figure. Selected N-factor contours are also superimposed to 
highlight the weaker harmonic disturbances. The computed frequency for the second-mode 
disturbance remains fairly co nstan t over the flare and falls within 2 % of the measured 266- 
kHz frequency at the base. The 266-kHz disturbance is thus identified as a second-mode 
disturbance. 

Amplification-rate spectra, -oti vs. F, are plotted in Fig. 3.32 for five stream wise 
locations from Re s = 2.6 x 10 6 to 4.5 x 10 6 . Electronic noise inhibits the detection of 
disturbance amplification at the upstream locations. Downstream, however, significant 
amplification is observed around the 15.6-kHz and 266-kHz frequency bands. These 
bands are attributed to the facility noise and the second-mode disturbance, respectively. 
Harmonics of the second-mode disturbance are amplified to a lesser degree at 133 kHz, 

531 kHz, and 797 kHz. 

Amplification-rate distributions, -0Ci vs. Re s , for the second-mode disturbance, the 
above harmonically-related disturbances, and the facility noise are plotted in Fig. 3.33. 

The facility noise experiences two regions of amplification which reach peak values around 
Re s = 3.3 x 10 6 and 3.9 x 10 6 . The downstream peak coincides with that of the second- 
mode disturbance. Small levels of sub-harmonic amplification are detected downstream of 
Re s = 3.5 x 10 6 . The first and second harmonics of the second-mode disturbance show 
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peak amplification at Re s = 4.1 x 10 6 and 4.3 x 10 6 , respectively. Transition onset was 
also detected in this vicinity, at Res = 4.2 x 10^, from surface-temperature measurements. 
The most-amplified of the second-mode harmonics is the second harmonic at 797 kHz 
whose peak-amplification rate is greater than twice that of the first harmonic and the sub 
harmonic. 

N-factor distributions, (N-No) vs. Res, are plotted for these frequencies in Fig. 

3.34. Facility noise grows slowly and reaches a plateau around Re s = 3.5 x 10 6 . Slightly 
downstream at Re s = 3.65 x 10 6 facility noise as well as the second-mode and sub- 
harmonic disturbances experience linear growth. In this region, the growth of the facility 
noise and the second-mode disturbance grow at the same rate. Sub-harmonic growth 
occurs at a much slower rate as indicated by the slope of the 133-kHz distribution. A small 
amount of growth of the second-mode harmonics is observed downstream of transition 
onset The first harmonic grows at a slower rate than the second harmonic, however. 

Fig. 3.35 shows the N-factor distributions for the second-mode disturbance 
measured at 266 kHz and for the LST-based computations of Balakumar and Malik 32 for 
the most-unstable disturbance at 220 kHz. The measured N-factor growth of the second- 
mode disturbance remains essentially linear, consistent with linear stability theory. The 
measured distribution was previously plotted in Fig. 3.34 as the difference, N - No, where. 
No is the N factor at the location where disturbance growth is initially detected. The value 
of No = 6.5 was determined from the LST computations. Electronic noise obscures the 
second mode from the most-upstream location up to Re s = 3.65 x 10 6 for this case. 
Downstream of this location, facility noise in the 0-50 kHz band does not inhibit detection 
of the 266-kHz disturbance. The 20% discrepancy between the measured and computed 
frequencies is due to the finite windward alignment angle of the test model. Beyond 
transition onset at Re s = 4.2 x 10 6 the slope of the measured N-factor distribution levels off 
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slightly. The deviation from linear growth may be accompanied by non-linear interactions 
with the second-mode harmonics. To determine if these non-linear phase coupling effects 
were present in the experiment the bispectrum was examined. 

The bispectrum magnitude, IB(Fi ,F2)1, is plotted as an image in Fig. 3.36 for the 
most-downstream measurement location, Re s = 4.5 x 10^. The Fi and F2 frequency axes 
are normalized by the second-mode frequency, F? n d Mode = 266 kHz, to facilitate 
identification of harmonically-related disturbances. An arbitrary logarithmic scale was used 
for the magnitude color map. Regions of deterministic phase coupling are indicated by 
localized peaks in the bispectrum magnitude. The diagonal line described by F2 = Fi is a 
line of symmetry, and the two regions on each side of this line contain identical 
information. Examining the triangular area along the vertical axis, phase coupling is 
observed for the following disturbances. The dark region where Fi = 266 kHz and F2 = 
266 kHz indicates that the first harmonic at F3 = (Fi + F2) = 532 kHz is phase coupled 
with the second-mode disturbance. Also, where Fi = 266 kHz and F2 = 531 kHz the 
bispectrum shows phase coupling of the second harmonic at F3 = 797 kHz with both the 
second mode and the first harmonic. A small degree of phase coupling is also observed 
along the line Fi + F2 = F2ndMode» where constituent disturbances are also frequency- 
coupled with the second mode. This line intersects the line of symmetry at the point Fi = 
F2 = 133 kHz, indicating some degree of phase coupling of the sub harmonic with the 
second mode. 

3.5.2 a = 2° Windward Case 

Power spectra are plotted in Fig. 3.37. The most-upstream spectrum looks very 
similar in shape to the electronic-noise spectrum. Two distinct peaks, however, develop 
downstream. The low-frequency peak is identified as facility noise in the 15.6-kHz band. 
The other disturbance is centered at 379 kHz. The normalized amplitude spectra are plotted 
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in Fig. 3.38. This plot shows more clearly growth in the region of 50-200 kHz that was 
obscured by electronic noise in the power-spectra plot The computed second-mode 
disturbance frequency under predicts the measured frequency by 1 1% as shown in the 
composite plot in Fig. 3.39. The discrepancy between these measured and computed 
frequencies is due both to the limited precision of the model alignment and to the limited 
accuracy of the X *= 28 approximation on which the computed estimate was based. 
Nevertheless, the 379-kHz disturbance is attributed to the second mode. 

At the most-upstream location. Res = 2.7 x 10 6 , significant amplification of the 
facility noise is observed in the amplification-rate spectra plotted in Fig. 3.40. At Re s = 4 x 
10 6 and 4.5 x 10 6 the second-mode, 379-kHz disturbance is also amplified. The 50-200 
kHz band is amplified to a lesser extent The amplification-rate distributions of individual 
disturbance modes are plotted in Fig 3.41. Facility noise is tracked by the 15.6-kHz and 
44.9-kHz modes. This noise is greatly amplified at the upstream measurement stations and 
moderately amplified downstream. The 143-kHz disturbance is representative of the 50- 
200 kHz band. This band experiences moderate levels of amplification over most of the 
flare which are slightly larger than those observed in this band for the zero-degree baseline 
case. Initial amplification of the second-mode, 379-kHz disturbance is detected just 
upstream of Re s = 4x 10 6 . A peak-amplification rate is reached and starts to decay near 
the end of the model The peak-amplification rate of the second-mode disturbance is 
significantly less than that observed in the baseline case, however. Slight amplification of 
the second-mode sub harmonic, indicated by the 189-kHz distribution, is observed near the 
base of the model. 

The growth of these modes is illustrated by the N-factor distributions plotted in Fig. 
3.42. The facility noise shows substantial growth from the initial location up to Re s = 3.5 
x 10 6 . Additional facility-noise growth is detected downstream of Re s = 3.75 x 10 6 . Near 
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this second region of growth, at Res = 4 x 10 6 , linear growth of the second-mode 
disturbance is detected. The low-frequency band represented by the 143-kHz mode shows 
linear growth over the entire flare region but at a much slower rate than the second-mode 
disturbance. Measured growth of the sub harmonic is not significant. 

The image of the bispectrum magnitude at Re$ = 4.5 x 10b, plotted in Fig. 3.43, 
shows only a slight degree of phase coupling. A band along the line Fi + ¥2 = F2nd Mode 
appears again, although not as prominent as in the baseline case. Since only moderate 
growth of the second-mode disturbance was measured and transition onset was not 
detected, strong non-linearities are not expected for this case. Because significant growth 
was not detected for any second-mode harmonics in the power spectra, localized phase 
coupling would also not be detected in the bispectrum. 

3.5.3 a = 4° Windward Case 

In addition to facility and electronic noise, the power spectra plotted in Fig. 3.44 for 
this case show evidence of intermittent mechanical resonance of the hot-wire sensor. 
Sporadic peaks appear at a few stream wise locations but are not characteristic of the 
physical flow. The hot-wire sensor used for this test case was near the end of its 
operational lifetime and, thus, prone to such resonance effects. However, detection of the 
second-mode disturbance is not inhibited, as shown by the normalized amplitude spectra 
image plotted with the computed second-mode frequency in Fig. 3.45. Significant growth 
is detected both in the facility-noise band and in the low-frequency band from 50-200 kHz. 
Amongst the sporadic resonance peaks, however, appears a smoothly growing disturbance 
which is characteristic of the physical flow. This smoothly-growing disturbance at 379 
kHz falls 3% below the computed second-mode frequency. The amplification-rate spectra 
plotted in Fig. 3.46 are degraded by the effects of the hot-wire resonance at several 
frequencies. There is significant amplification in the facility-noise and low-frequency 
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bands. The second-mode disturbance at 379 kHz also shows significant amplification 
downstream. 

Amplification-rate distributions for frequencies not corrupted by hot-wire resonance 
are shown in Fig. 3.47. Facility noise is indicated by the 15.6-kHz and 44.9-kHz 
distributions, and the low-frequency band from 50-200 kHz is represented by the 105-kHz 
distribution. The peak-amplification rate of the low-frequency band has increased relative 
to the 2-degree windward case. The facility noise and the 105-kHz disturbance all reach a 
maximum amplification at Res = 3.7 x 10 6 . At this same location, amplification of the 
second-mode disturbance is initially detected, as indicated by the 379-kHz distribution. 

The peak-amplification of this second-mode disturbance has further decreased relative to 
the 2-degree windward case. The sub harmonic, at 189 kHz, is amplified by a small 
am ount downstream of peak in the second-mode amplification. The N-factor distributions 
plotted in Fig. 3.48 show significant growth of the facility noise as well as the following 
features. Growth of the low-frequency disturbance at 105 kHz is detected downstream of 
Re s = 3.5 x 10 6 . Growth of the second mode is detected downstream at Re s = 3.8 x 10 6 
but occurs at a faster rate. Measured growth of the sub harmonic is not significant. 

3.5.4 a = 2° Leeward Case 

In contrast to the previous cases, the upstream power spectrum for this case, shown 
in Fig. 3.49, bears little resemblance to the electronic-noise spectrum. Significant flow- 
disturbance growth has already occurred prior to the most-upstream measurement location. 
Initially, flow disturbances are primarily confined to the frequency range 0-200 kHz. A 
peak-power disturbance is initially detected at 141 kHz, but subsequently disperses into a 
more broad-band disturbance downstream. Spectral broadening is also illustrated by the 
normalized amplitude spectra in Fig. 3.50. Here, measured growth is initially restricted to 
the 0-200 kHz band. Downstream, however, growth is observed over a wider band up to 
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600 kHz. Near the base the rate of growth decreases, indicating saturation takes place. A 
composite plot of the computed second-mode frequency and the amplitude-spectra data is 
shown in Fig. 3.51. The computed frequency is about 45 kHz at the base. This mode is 
not clearly identified from the spectra due to the broad-band character of the flow 
disturbances. In addition, the second mode is obscured because its frequency falls within 
the 0-50 kHz facility-noise band. However, the initial peak-power disturbance falls within 
the 50-200 kHz band and, thus, is not associated with the second mode. 

Amplification-rate spectra, plotted in Fig. 3.52, confirm that the most-unstable 
disturbances fall below 200 kHz at the initial upstream location. At downstream locations, 
the higher-frequency disturbances from 200-600 kHz experience a moderate increase in 
amplification indicating that spectral broadening is taking place. Amplification rates decay 
to negligible values at the final streamwise location as saturation occurs. Facility noise is 
again represented by the 15.6-kHz and 44.9-kHz disturbances for the amplification-rate 
and N-factor distributions presented in the following. Also, the 50-200 kHz band is 
represented by both the 60.5-kHz and 141-kHz distributions. The development of the 449- 
kHz disturbance is characteristic of the high-frequency band from 200-600 kHz. 

Amplification-rate distributions of the above representative modes are plotted in 
Fig. 3.53. Amplification of the facility-noise and low-frequency modes remains significant 
up to Re s = 4 x 10 6 . Amplification of the 141-kHz disturbance decays quickly, however, 
and is even attenuated in the vicinity of Res = 3.6 x 10^. Here, attenuation is not attributed 
to linear damping but rather to non-linear disturbance interaction. This non-linear 
interaction is demonstrated by amplification of the 449 -kHz, representative disturbance. As 
the 141-kHz dominant mode is attenuated, the 449-kHz mode is amplified. This reveals an 
amplification shift from the 50-200 kHz band to the 200-600 kHz band, an indication of 
non-linear frequency leakage. 
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N-factor distributions, plotted in Fig. 3.54, illustrate linear growth and subsequent 
non-linear saturation of the representative modes. For the 141-kHz, 60.5-kHz, and 
facility-noise modes the initial growth is approximately linear. The 141-kHz disturbance, 
however, saturates at the transition-onset location, Re s = 3.33 x 10 6 . Downstream, at Res 
= 4 x 10 6 , the low-frequency and facility-noise modes saturate. These saturated modes 
stimulate growth of higher-frequency disturbances, represented by the 449-kHz mode. 

The 449-kHz distribution shows steady continuous growth downstream of the transition- 
onset location. 

The coincidence of the saturation point of the 141-kHz disturbance and the 
transition-onset location suggests that this mode is the dominant instability mechanism. 
Since this frequency falls much higher than the computed second-mode frequency and the 
approximate range of first-mode disturbances, the dominant mode must be attributed to a 
different physical mechanism. Another possibility for the dominant mechanism on the 
leeward ray is the cross-flow instability. This instability may arise when the velocity 
component normal to the edge velocity reaches a maximum value within the boundary 
layer.38 The origin of this instability for the flared-cone model is illustrated by the surface 
streamlines and velocity vectors shown in Fig. 3.55 obtained from Navier-Stokes 
computations for the a = 2° configuration. The angle of attack, a, shown in the X-Z plane 
is greatly exaggerated for clarity. The diametrically opposing windward and leeward rays 
are indicated in the figure. Surface streamlines were actually evaluated on the K = 2 
surface from the numerical mesh since the K = 1 surface is a no-slip boundary. These 
laminar computations predict a small separated region along the leeward ray, just 
downstream of the tangent Separation was not observed in the a = 2° leeward case, 
however. Velocity vectors at an arbitrary location along the lateral meridian in the X-Y 
plane, at X = 5", which corresponds to Re s = 1.17 x 10 6 and is referred to as point A, are 
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shown in a detailed view in the figure. Looking directly upstream from the point of view 
indicated the model surface is on the right and the farfield boundary is on the left The 
undisturbed freestream velocity indicated by the circles is in the outward direction, normal 
to the plane of the paper. Across the shock the velocity is deflected away from the surface. 
Within the boundary layer the velocity decreases and changes direction to follow surface 
streamlines around the model. The peak in the cross-flow velocity is readily apparent in 
this vector plot Disturbances originating from this instability propagate along local 
streamlines from the windward ray towards the leeward ray. Thus, the dominant instability 
on the leeward ray was attributed to a cross-flow type of instability, however, no direct 
measurement of the frequency of such a cross-flow disturbance was obtained. 

Fig. 3.56 shows the image of the bispectrum magnitude at transition onset, Re s = 
3.33 x 10 6 . The frequency axes are normalized by the frequency of the cross-flow 
disturbance, FcrossFlow = 141 kHz. Although the cross-flow disturbance is at its 
maximum amplitude, the spectrum is fairly broad band at this location. The bispectrum 
shows a very small amount of broad-band phase coupling in the low-frequency band. The 
absence of strong phase coupling in localized bands indicates that disturbances are random 
and chaotic. However, a slight degree of localized phase coupling is observed along the 
line Fi + F2 = Fcross Flow- This localized band is much less distinguished, however, than 
the si milar features observed in the a = 0° baseline and a = 2° windward cases. 

3.5.5 a = 4° Leeward Case 

Separation was detected from surface-temperature measurements for this case. The 
following measurements are therefore presented to illustrate the development of turbulent- 
flow disturbances through the bypass-transition process. Power spectra plotted in Fig. 
3.57 show the last stages of spectral broadening taking place at the most-upstream 
locations. Distinct peaks are detected at 15.6 kHz and 82 kHz which fall within the 0-50 
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kHz facility-noise and the 50-200 kHz low-frequency bands, respectively. The peak- 
power disturbance is at 141 kHz, which may have saturated upstream of the initial 
measurement location. These distinct peaks coalesce downstream as the spectrum 
broadens. The image of the normalized amplitude spectra, shown in Fig. 3.58, does not 
clearly identify the second mode. The computed second-mode frequency is around 28 
kHz, which falls within the facility-noise band. Growth is measured in the facility-noise 
and low-frequency bands as well as in the 200-600 kHz band to a lesser extent. The 
saturated 141-kHz disturbance does not grow any further in the stream wise direction and 
appears as a narrow-band horizontal stripe in this image. Amplification-rate spectra are 
plotted in Fig. 3.59. At Res = 2.6 x 10 6 amplification is observed in the facility-noise band 
and in the high-frequency band from 200-600 kHz. The 141-kHz disturbance is attenuated 
at this location as a result of growth saturation. At downstream stations amplification 
becomes negligible. As in the a = 2° leeward case, the dominant disturbance is the 141- 
kHz mode which is also attributed to a cross-flow type of instability mechanism. 

Amplification-rate distributions of discrete modes are plotted in Fig. 3.60. The 
facility noise is again represented by the 15.6-kHz mode. The estimated second-mode 
disturbance frequency, indicated by the 27.3-kHz distribution, falls within the facility-noise 
band. Disturbances within the 50-200 kHz band are indicated by the 82-kHz and 141-kHz 
distributions. In addition, the 200-600 kHz band is represented by the 449-kHz mode. 
With the exception of the 141-kHz disturbance, the largest amplification occurs at the most- 
upstream measurement station. The saturated 141-kHz disturbance is initially attenuated. 
Downstream of Re s = 3 x 10 6 flow disturbances experience small levels of amplification. 
The net growth is illustrated in the N-factor distributions shown in Fig. 3.61. The 82-kHz 
and 141-kHz modes show negligible growth as both remain fairly constant over the flare 
confirming that saturation has taken place prior to the initial measurement station. 
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Saturation of the high-frequency 449-kHz mode is observed slightly downstream at Re s 
3 x 10 6 . The facility noise continues to grow moderately and finally saturates near the 

base. 

3.6 Summary of Angle-of-Attack Effects 

The above stability and transition data as well as the dominant-instability 
mechanisms which lead to transition are summarized below. The presentation of results is 
then concluded with a discussion of the angle-of-attack effects on the boundary-layer 
stability diagram. 

3.6.1 Stability & Transition Data 

The stability and transition data are listed in Table 3.1. The local Reynolds number 
is listed for the locations where maximum amplification, transition onset, and growth 
saturation were observed. In Fig. 3.62. local Reynolds numbers normalized by the 
transition-onset Reynolds number at a = 0° are plotted against the ratio of the angle of 
attack to the cone half angle, (Xcone = 5°. Thus, on the horizontal axis, 2 and 4 are 
indicated by 0.4 and 0.8, respectively. A positive value refers to the windward ray, and a 
negative value represents the leeward ray. The Reynolds numbers at the first and last hot- 
wire measurement stations are indicated by the dashed and dashed-dotted lines, 
respectively. Plot symbols that fall outside the measurement region indicate approximate 
locations. For the zero-degree baseline case transition onset is preceded by the maximum 
amplification of the dominant disturbance, labeled Maximum l-oql on the plot. Growth 
saturation, indicated by the Maximum (N-Nq) symbol, was not observed within the 
measurement region, thus, the exact location is unknown. 

Transition onset is delayed for both windward cases. Transition and saturation 
were not observed in these cases, but the maximum-amplification point demonstrates the 
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asymptotic stabilizing trend. The stabilizing effect is somewhat reduced for the 4-degree 
windward case. Transition onset shifts upstream for the leeward cases, relative to the 
baseline case. For the 2-degree leeward case transition onset is also accompanied by 
saturation. This coincidence suggests that the combined effects of adverse-pressure 
gradient and angle of attack accelerate the process of non-linear breakdown for the 
dominant instability on the leeward ray. Transition onset for the 4-degree leeward case was 
deter mine d by the measured separation point Here, separation provides the bypass 
mechanism for instantaneous transition to turbulence. However, maximum amplification 
and saturation of the dominant disturbance occurred upstream of the measured separation 
location. Thus, the natural-transition process had progressed significantly prior to the 
observed separation point 
3.6.2 Dominant-Instability Mechanism 

Present computations show that angle of attack has an asymptotic effect on the 
second-mode disturbance frequency. This effect is graphically illustrated in Fig. 3.63 
where the frequency normalized by the a = 0° value are plotted versus the normalized angle 
of attack. However, the second-mode disturbance is not the dominant-instability 
mechanism in all test cases. Dominant instabilities observed in the leeward cases were 
attributed to a cross-flow type of instability and are indicated as such on the plot The 
dominant disturbance frequency was 141 kHz for both the leeward cases. The facility- 
noise floor delineates the 0-50 kHz band in which flow disturbances were not uniquely 
identified in the experiment The a = ±4° cases followed similar trends as the a = ±2° 
cases, thus stability diagrams are presented in the following section only for the a = 0° and 
±2° cases. 


48 



3.6.3 Stability Diagrams 

A stability diagram consists of amplification-rale contours which identify the 
regions of disturbance frequencies and spatial locations that are stable and unstable. In the 
stable region flow disturbances are attenuated and ~cq < 0. The unstable region is the 
region where -oti > 0 and disturbances are amplified. The neutral-stability curve separates 
these two regions and is equivalent to the amplification-rate contour where -cq = 0. The 
present stability diagrams were constructed by plotting the measured amplification-rate data 
as an image with F vs. Re s . All regions where -oq^Ointhe measured stability diagram do 
not necessarily correspond to stable regions of the theoretical stability diagram, however. 
Such discrepancies between measured and theoretical stability diagrams are nonetheless 
expected for the following reasons. Amplification is not detectable until the flow- 
disturbance signal exceeds the noise level in the hot-wire anemometry system. Thus, initial 
regions of amplification were not experimentally detectable. Additionally, in some cases, 
downstream regions of apparent neutral stability were attributed to non-linear processes 
such as saturation and mean-flow distortion, Le., transition onset. In spite of these 
limitations, however, the measured stability diagrams are a useful tool in graphically 
illustrating angle-of-attack trends. 

The stability diagram for the a = 0° baseline case in the frequency range 0-500 kHz 
is shown in Fig. 3.64. Most-unstable regions are indicated by the areas in the image where 
-oti has largest values which are denoted by the color map levels on the right of the plot. 
The facility-noise and second-mode disturbances in the 0-50 kHz and 266-kHz bands, 
respectively, are the most-unstable disturbances. Narrow-band harmonics of the second 
mode at 532 kHz and 797 kHz, not shown on this diagram, were also amplified. 
Amplification in the 50-200 kHz band is also detected. Fig. 3.65 shows the stability 
diagram for the a = 2° windward case. The 379-kHz second-mode disturbance is clearly 
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seen. Moderate amplification is again apparent in the 50-200 kHz band. Amplification of 
the facility noise is shifted upstream because the physical location of the model within the 
test nozzle varies with angle of attack. Broad-band disturbances are apparent in the stability 
diagram of the a = 2° leeward case, shown in Fig. 3.66. Non-linear saturation of the 141- 
kHz dominant disturbance attributed to a cross-flow type of instability is indicated by the 
region where -0Ci= 0 at Re s = 3.33 x 10 6 . The maximum-amplification point of this 
disturbance occurred upstream of the initial measurement station. 

A combined stability diagram is shown in Fig. 3.67 for the a = +2°, 0°, and -2° 
cases. Amplification-rate contours are plotted in the range -otp 4-9 x 10"^. These contour 
levels were chosen to highlight the maximum-amplification point of the dominant 
instabilities. Although the neutral-stability curve was not accurately determined, the 
maximum-amplification point is also a significant feature of the stability diagrams. This 
composite plot clearly demonstrates the angle-of-attack effects on both the frequency and 
the region of amplification for the dominant instabilities, with respect to the ct = 0° baseline 
case. For the windward case the frequency increases by 43%, and the Reynolds number at 
maximum amplification is shifted downstream by 9.4%. In contrast, the maximum- 
amplification Reynolds number is shifted upstream by 40% for the leeward case, and the 
frequency decreases 47%. Because the angle-of-attack effects are asymptotic, the a = ± 4° 
cases are essentially coincident with the a = ± 2° cases and are omitted for clarity. 
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4 Concluding Remarks 


The first investigation of angle-of-attack effects on the stability of the hypersonic 
boundary layer on a conical model in a low-disturbance facility has been conducted. The 
experiments were performed in the Nozzle-Test-Chamber Facility at NASA Langley 
Research Center. The 5° half-angle flared-cone model remained within the low-noise 
uniform-freestream envelope of the Mach-6 Axisymmetric Quiet Nozzle for the 0°, 2°, and 
4° angles of attack investigated. The freestream Mach number for the present tests was M« 
= 5.91 and the freestream Reynolds number per unit length was Re/1 = 2.82 x 10 6 ft -1 . 
Laminar-flow solutions were obtained from Navier-Stokes computations to model the mean 
flow over the flaied-cone geometry. The adverse-pressure gradient was documented with 
surface static-pressure measurements. Transition onset was detected with surface- 
temperature measurements. The state of the boundary layer was verified by schlieren flow 
visualization. Constant-voltage anemometry diagnostics were applied to obtain both mean- 
flow and stability data from a series of point measurements with a single-component, hot- 
wire probe. Amplification rates and growth factors of normal-mode disturbances within 
the boundary layer were measured. Higher-order spectral analysis was applied to detect 
non-linear disturbance interactions. 

4.1 Summary of Results 

Significant findings of the present work are summarized as follows. First, the 
effects of the adverse-pressure gradient on the flared-cone model at a zero-degree angle of 
attack are discussed. Then, the effects of angle of attack are summarized. Finally, the 
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effects of facility noise on the stability and transition measurements in the low-disturbance 
wind tunnel are addressed. 

4.1.1 Effect of Adverse-Pressure Gradient 

Transition occurred at a Reynolds number of 4.2 x 10 6 on the flared-cone model 
for the zero-degree angle-of-attack case. In contrast, on straight-walled cones transition 
Reynolds numbers may be in excess of 10 x 10 6 . The natural-transition process was 
accelerated by the adverse-pressure gradient over the flare. The adverse-pressure gradient 
increased amplification rates, relative to a zero-pressure gradient cone, for both first and 
second-mode disturbances. 32 The adverse-pressure gradient also retarded boundary-layer 
growth over the flare. The frequency of the second-mode disturbance decreases with the 
increasing boundary-layer thickness on a straight cone. 16 However, the slowly decreasing 
boundary-layer thickness over the flare resulted in no significant change in the measured 
second-mode frequency in the present tests. The measured location of maximum- 
disturbance energy in the boundary layer corresponded to the computed generalized- 
inflection point 

The measured rate of growth of the second-mode disturbance compared well with 
the LST computations of Balakumar and Malik 32 in regions where the signal-to-noise ratio 
of the second mode was greater than 1. A discrepancy between the predicted and measured 
second-mode frequencies was attributed to the finite alignment angle of the model in the 
experiment Transition onset is quickly followed by non-linear breakdown of the dominant 
second-mode disturbance. Non-linear phase coupling between the dominant second-mode 
disturbance and harmonically related disturbances was observed. The peak-amplification 
rate of the second harmonic exceeded twice that of the sub harmonic. 
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4.1.2 Angle-of-Attack Results 

On the windward ray the following observations were made. The boundary layer 
remained laminar over the length of the model. Thus, the transition Reynolds number 
exceeded the length Reynolds number of 4.7 x 10 6 . The frequency of the dominant 
second-mode instability increased relative to the zero-degree case, scaling with the thinner 
boundary-layer thickness. Also, growth and amplification of the second-mode instability 
were detected at larger local Reynolds numbers. The peak-amplification rate of the second- 
mode instability decreased significantly with increasing angle of attack. Amplification of 
the low-frequency band, associated with the first-mode instability, slightly increased with 
increasing the angle of attack. But, these disturbances did not grow sufficiently to induce a 
mean-flow dis tortion. Increasing the angle of attack from 2° to 4 did not measurably 
change the frequency of the second-mode instability. 

Along the leeward ray the following trends were noted. Transition occurred at a 
local Reynolds number of 3.33 x 10 6 for the 2-degree case. Transition was induced by 
separation at 2.75 x 10 6 for the 4-degree case, however, significant normal-mode 
disturbance growth had occurred upstream of the observed separation point. The estimated 
second-mode frequency decreased to well below the frequency of the dominant-instability 
mechanisms. The dominant disturbance was attributed to a cross-flow type of instability 
originating on the side of the model. The frequency of the dominant instability remained 
constant as the angle of attack increased from 2° to 4°. The peak-amplification rate occurred 
at lower local Reynolds numbers than in the zero-degree case. Non-linear growth 
saturation and spectral broadening were observed in the transitional and turbulent regions. 
4.1.3 Effect of Facility Noise 

Within the low-noise test volume, disturbance growth was measured in the 0-50 
kHz facility-noise band. This measured growth corresponded to rising levels of noise in 
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the stream wise direction, however, as opposed to normal-mode instabilities. The 
occurrence of growth in the facility-noise band did not correlate spatially with the location 
of transition onset in any test case. In addition, angle of attack did not affect the overall 
measured growth of this band. The shift in location of measured facility-noise growth was 
attributed to the change in physical position of the model at angle of attack within the test 
nozzle. Furthermore, the low-frequency facility-noise band did not obscure detection of 
the high-frequency dominant instability mechanisms. The present higher-order spectral 
analysis did not reveal any deterministic phase coupling with the dominant instabilities and 
the facility-noise band. 

4.2 Conclusions 

• Second-mode instabilities dominate the natural-transition process on the flared-cone 
model at a zero-degree angle of attack. 

• The second-mode frequency is highly sensitive to a small angle of attack. However, the 
change in the second-mode frequency with increasing angle of attack is asymptotic. 

• The second mode was stabilized along the windward ray, but the adverse-pressure 
gradient destabilized the first mode. 

• Along the leeward ray the estimated second-mode frequency decreased well below the 
dominant-instability mechanism band associated with transition. 

• Transition on the leeward ray may be due in part to disturbances originating from a cross- 

flow type of instability. 

• The low-noise disturbance field in the test facility did not prematurely induce transition. 
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4.3 Future Work 


Recommendations for future work are discussed below in the general areas of 
stability experiments, stability computations, and flow visualization. 

4.3.1 Stability Experiments 

Because the present stability measurements were obtained with a single-component 
hot-wire probe, no information about the orientation of flow disturbances was available. 
The dominant second-mode disturbances are most unstable when propagating along the 
mean-flow direction, but first-mode disturbances are most unstable when propagating at an 
oblique angle. The application of multiple-probe correlation techniques recently 
demonstrated by Kimmel et al 39 is recommended for providing the detailed information of 
the disturbance-field orientation. Additionally, the application of such techniques is 
recommended to quantify the influence of freestream disturbances on boundary-layer 
stability and transition with more certainty. The need for understanding this influence is 
illustrated by the following. 

For a conical model at angle of attack the second-mode frequency spans a wide 
range of values from the windward ray to the leeward ray. Also, first-mode instabilities 
and cross-flow insta bilitie s are present. The freestream-disturbance field may interact with 
any of these mechanisms. To assess the influence of this interaction on transition 
measurements, Reed and Haynes 40 presented transition correlations for regions of cross 
flow on a cone at an angle of attack from experiments conducted in both “quiet” 9 and 
“noisy” 10 wind-tunnel facilities. Recent transition measurements 8 compared well with the 
“quiet” correlation in regions of larger cross flow velocities, although these tests were 
conducted in a “non-quiet” shock-tunnel facility. The explanation for this unexpected 
agreement was given that the facility-disturbance field did not interact with the dominant 
cross-flow instabilities on the cone at angle of attack. However, interaction of the facility- 
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disturbance field with other instability mechanisms present in the experiment in regions 
where cross-flow velocities are small was not assessed. With multiple-probe 
measurements and the application of spatial and temporal correlation techniques, the 
receptivity of the model boundary layer to the freestream-disturbance field can be quantified 
and ultimately better understood. 

4.3.2 Stability Computations 

Stability computations for the angle-of-attack configuration were not available to 
identify specifically the first and second modes. Thus, angle-of-attack trends were based 
on observations of the representative modes which were typical of normal-mode 
disturbances in each respective frequency band. These angle-of-attack trends are unique to 
the flared-cone model, however. For example, the maximum amplification decrease for the 
second mode and increase for the first mode on the windward ray were not observed in 
previous straight-cone angle-of-attack experiments.^ Therefore, to better understand the 
combined effects of angle of attack and adverse-pressure gradient on instability growth and 
the natural-transition process, a computational stability study is recommended. The 
database of measurements from the present stability experiments conducted in a low- 
disturbance wind tunnel are suited for validation of stability theory computational tools. 
Also, stability computations are useful in the interpretation of stability measurements. 

Application of the LST-based e N method to the zero-degree configuration has been 
previously demonstrated by Balakumar and Malik.32 Computational effort for an angle-of- 
attack study can be reduced considerably by having a guide as to the range of frequencies 
and locations over which disturbances are unstable.^ To that end, the present measured 
stability diagr ams can be used as such a guide in judiciously choosing the frequencies and 
locations over which stability computations are performed. 
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While the eF method is widely used for stability computations, significant 
stream wise and circ umf erential variation in the mean-flow over the flared-cone model at an 
angle of attack may violate the parallel-flow approximation on which this LST method is 
based. However, emerging methods based on the Parabolized Stability Equations (PSE) 
have shown promise for modelling the stability of the non-parallel flow over a cone at angle 
of attack. 4 ! jn addition, the PSE-based method may be used to model certain non-linear 
disturbance interactions such as those observed in the present experiments 4 ^ Based on 
these considerations, application of the PSE-based method is also recommended for the 
computational stab ilit y analysis of the flared-cone geometry at angle of attack. In 
performing any of the above stability computations special consideration should be given to 
the formulation of a definition of the disturbance N-factor in a three-dimensional non- 
parallel flow. 

Mean-flow solutions for stability computations may be obtained with one of several 
computational tools such as the CFL3D code used in the present work. Discrepancies 
observed in the angle-of- attack cases between the computed adiabatic-wall temperature and 
the measured surface temperature, however, indicated that thermal conduction effects were 
present in the experiment. Stability computations based on such an adiabatic mean-flow 
solution would not account for these conduction effects. Application of the measured 
steady-state temperature as a boundary condition in the mean-flow computations, however, 
would properly simulate the thermal-equihbrium state of the test model. Thus, obtaining 
measurements of the circumferential temperature distribution is recommended to establish 
this boundary condition. Unfortunately, the present measurements were obtained only in 
the symmetry plane. To obtain the circumferential measurements the test model must be 
rotated about its axis. Rotating the model was deemed damaging to delicate thermocouple 
wires and pressure tubing extending through the model sting. The thermocouple wires are 
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easily fatigued, and the pressure-tubing connections may leak when over stressed. 
Therefore, in order to preserve the integrity of the internal model instrumentation the 
circumferential temperature distribution was not measured in the present angle-of-attack 
study. 

4.3.3 Flow Visualization 

A compilation of earlier studies in which hypersonic boundary-layer instability 
waves have been observed in wind tunnel tests was presented by Smith. 4 ^ These waves 
have been characterized as second-mode disturbances whose wavelength is approximately 
twice the boundary-layer thickness. Visualization methods applied in these studies include 
schlieren and shadowgraph photography. These methods integrate density changes along a 
collimated light path, thus, no three-dimensional visualization of the instability waves is 
currently available. Detailed computations of Pruett and Zang 44 using Direct Numerical 
Simulation predicted that the instability wave is a three-dimensional disturbance with 
significant circumferential variation. A three-dimensional visualization method is needed 
for validation of this numerical result 

One such method suited for three-dimensional visualization of high-speed flows in 
wind tunnels is the laser-holographic focusing-schlieren technique. An improved system 
was proposed by Weinstein 45 and first experimentally demonstrated by Doggett and 
Chokani 46 in the Supersonic Wind Tunnel Facility at North Carolina State University. The 
laser-holographic focusing-schlieren method can record an entire flow field in a single 
hologram during the short-duration pulse of the laser light source. The hologram can then 
be used to reconstruct detailed images of the flow in the test section. An original goal of 
the present work was to install the laser-holographic focusing-schlieren system developed 
at NCSU in the Nozzle-Test-Chamber Facility at NASA Langley. Unfortunately, due to 
misalignment during shipment the pulsed laser could not be restored to its peak 
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performance during the limited amount of time allotted for this temporary installation. 
Future application of the laser-holographic focusing schlieren technique is recommended to 
validate the predicted three-dimensional structure of hypersonic boundary-layer instability 
waves. 
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6 Tables 


Table 2. 1. Hot-wire probe dimensions 


No. dp dg dl 


1 



0.0185 

185 

2 



0.0215 

215 

3 

0.005 

0.025 

0.021 

210 


Notes: All dimensions are in inches except for the dimensionless length-to- 

diameter ratio. 


Table 2.2. Test cases & measurements 


a 

Mean-Flow 

Surface 

Pressure 

Mean-Flow 

Surface 

Temperature 

Mean-Flow 

Boundary 

Laver 

Boundary- 

Layer 

Stability 

Schlieren 

Flow 

Visualization 

44° 


X 

X 

X 

X 

+2° 

X 

X 

X 

X 

X 

CP 

X 

X 

X 

X 

X 

-2° 

X 

X 

X 

X 

X 

4° 


X 

X 

X 

. - ■ 1 

X 


Notes: x indicates measurement was obtained. 
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Table 3.1. Stability & transition data 


Reynolds Numbers at: 


Frequency (kHz) ofT 


Maximum Transition 
Amplification Onset 


Maximum 

Growth 


Measured Computed 

Dominant Second-Mode 

Instability Instability _ 


44 ° 

+ 2 ° 

CP 

- 2 ° 

4° 


4.1 

4.31 

3.94 

< 2.6 

< 2.6 


>4.5 

>4.5 

4.2 

3.33 

2.6 


>4.5 

>4.5 

>4.5 

3.33 

< 2.6 


379 

379 

266 

141 

141 


390 

337 

262 

45 

28 


Notes* A Dositive angle of attack denotes the windward ray. Local Reynolds numbers are 
based on s^ace^uc kngthand freestream conditions. Computed second-mode frequency 
is based on the X « 2 8 approximation and Navier-Stokes computations. 


65 




7 Figures 




layer 

Fig. 1.2. Low-disturbance wind-tunnel design 
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Fig. 2.2. Uniform Mach-6 Envelope and Quiet-Flow Regions in the NTC Facility. 
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Fig. 2.3. Geometry of Flared-Cone Test Model. 
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Fig. 2.5. Photograph of the Traverse Mechanism and Hot-Wire Probe Support in the 
NTC Facility with the Flared-Cone ModeL 



Fig. 2.6. Photograph of the Short-Strut Mount and the New Mounting Blocks 
Designed for the Present Angle-of- Attack Study. 
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Fig. 2.7. Photograph of Rared-Cone Model Installation with New Mounting Blocks 
in NTC Facility. 
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Hidden-line plot of the numerical grid for the flared-cone model geometry 
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Fig. 3.6. Measured and computed surface pressure for a = 0° baseline case with 
computations for a = 0.2°. 
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Fig. 3.7. Measured and computed surface pressure for a = 2° windward and 
leeward cases. 
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Fig. 3.8. Measured and computed surface temperature for a = 0° baseline case with 
computations for a = 0.2°. 
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Fig. 3.9. Measured and computed surface temperature for a = 2° windward and 
leeward cases. 
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Fig. 3.10. Measured and computed surface temperature for a = 4° windward and 
leeward cases. 
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Fig. 3. 1 1 . Schlieren flow visualization for a = 0° baseline case. 



Fig. 3.12. Schlieren flow visualization for a = 2° windward case. 
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Fig. 3.13. Schlieren flow visualization for a = 4° windward case. 
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Fig. 3. 14. Schlieren flow visualization for a = 2° leeward case. 
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Fig. 3.15. CV A mean-output profiles for a = 0° baseline case. 



Fig. 3.16. Measured and computed boundary-layer thickness for a = 0° baseline case 
with computations for a = 0.2°. 
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Fig. 3.18. CVA mean-output profiles for 










Fig. 3.23. Measured and computed generalized-inflection point for a - 0° baseline 
case with computations for a = 0.2°. 
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Fig. 3.29. Measured power spectra for a = 0° baseline case. 
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Fig. 3.34. 


Fig. 3.35. 
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computations of Balakumar & Malik. 32 
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Fig. 3.37. Measured power spectra for a = 2° windward case. 



Fig. 3.38. Normalized amplitude spectra for a = 2° windward case. 
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Fig. 3.44. Measured power spectra for a - 4° windward case. 
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Fig. 3.45. Computed second-mode disturbance frequency with measured spectra for 
a = 4° windward case. 


















15 



Fig. 3.52. Amplification-rate spectra for a = 2° leeward case. 
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Fig. 3.54. 
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N-factor distributions for a = 2° leeward case. 
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Fig. 3.55. Computed surface streamlines and velocity vectors for a - 2° 
configuration. 
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Bispectrum magnitude for a = 2° leeward case at Re s = 3.3 x 10 6 . 







Fig. 3.57. Measured power spectra for a = 4° leeward case. 
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Fig. 3.58. 
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spectra for a = 4° leeward case. 
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Fig. 3.6 1 . N-factor distributions for a = 4° leeward case. 
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Fig. 3.62. Maximum-amplification, transition-onset, and growth-saturation Reynolds 
numbers versus angle of attack. 
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Fig. 3.63. Computed second-mode disturbance frequencies at the model base versus 
angle of attack. 
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Fig. 3.64. Stability diagram based on measured amplification-rate spectra for a « 0° 
baseline case. 
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8 Appendices 


8.1 Angle-of- Attack Mounting Blocks 

A set of four mounting blocks were designed for the present angle-of-attack 
experiments conducted in the NTC. Blocks for both angle-of-attack and sideslip 
configurations were constructed to obtain schlieren and hot-wire measurements as 
discussed in section 2.1.7 of the Approach chapter. Fig. 8.1 below shows the position of 
the model inside the nozzle with both the sideslip and the angle-of-attack mounts for the 
maximum 4-degree angle investigated. This figure shows the model located within the 
low-noise uniform freestream of the “quiet” nozzle. Three-view drawings for each block 
are next shown in Figs. 8.2-8.5. 
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Fig. 8.1. Top and side views of flared-cone model installation in the NTC Facility. 
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Fig. 8.3. Three-view drawing of 2-degree angle-of-attack block, 
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Fig. 8.4. Three-view drawing of 4-degree angle 
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8.2 Hot-Wire Coordinate Transformation 


This appendix describes the coordinate transformation applied to the reference 
contact switch coordinates as explained in section 2.4.1 of the Approach chapter. The 
location of the contact switch was measured with respect to the nozzle coordinate system. 
The pivot point of the angle-of-attack mounting blocks was taken as the origin. The X 
coordinate was in the downstream direction along the nozzle center line, and Y was 
outward from the center line in the plane of the hot-wire measurements. 

First, the location of the hot wire was determined from the contact-switch 
coordinates and the probe-tip dimensions. The probe-tip dimensions, dp and dq, and the 
angle of the probe arm, 0, which are shown in Fig. 8.6, determine this coordinate shift as 
rshift = -(dq cos 0 + dp sin 0) i + (dq sin 0 + dp cos 0) j (8.1) 

where bold typeface denotes a vector quantity, and i and j represent unit vectors in 
the X and Y directions, respectively. The hot-wire location in nozzle coordinates was then 
referenced to the model coordinate system. This procedure is explained below, but first the 
model coordinate system and geometry are defined. 

The origin of the model coordinates is located at the sharp tip as shown in Fig. 8.7 
where x is the axial distance from the tip and y is radially outward in the measurement 
plane. The fore body of the model is a straight-walled cone of half angle Qcone • The 
profile of the aft section is a circular arc of radius rflare • This flared surface is tangent to the 
cone at xtangent, y tangent- The center of the circular-arc flare is located at x/iare, y flare as 
shown in the figure and evaluated by 

Xflare ~ Xtangent ~ * flare SW Qcone 

yflare = Xtangent tQJl Qcone + Xflare COS Qcone (8.2) 

The arc length along the flare surface is given by the following 


115 



where the reference 


COS Ocone 


angle 0, shown in Fig. 8.7, is given by 

0 = cos'i * ~ x farA 

T l r flare i 


( 8 . 3 ) 


(8.4) 


The local slope of the flared surface is 

dy _ ( X’Xflare ) 

dx (y - y flare) (8.5) 

At a zero-degree angle of attack the cone coordinates are related to the nozzle 

coordinates by the following 

X = Xexit + X - Xbase + ( Xbase ~ %exit) 

Y = y (8.6) 

where Xexit is the nozzle coordinate at the nozzle-exit plane and xbase is the cone coordinate 
at the model base. The axial distance between the model base and the nozzle-exit plane, 
(xbase - Xexit), was measured with the model mounted in the test section as illustrated in 


Fig. 8.8. 

The location of the model at an angle of attack can be described by the following. 
Let Xa, Ya represent the nozzle coordinates of a point A on the model surface, xa, yA» at a 
zero-degree angle of attack. In polar coordinates this point corresponds to Ra, 

When placed at an angle of attack, ex, point A is rotated to location B in the nozzle, Xb, Yb, 
which is specified by the polar nozzle coordinates Ra, The rotated coordinates 


are thus given by 

Xb = Ra cos [Qa + a ) 

Yb = Ra sin ( Qa + a) (8.7) 

The coordinate transformation in Eq. 8.6 requires a value for {xbase - XexitX 
However, when the model was mounted at an angle of attack the position of the model in 
the nozzle was determined by measuring (Xbase ■ Xexit)* Thus, an iterative procedure was 
applied to determine a value for (xbase - Xexit) such that when xbase was transformed and 
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rotated, via Eqs. 8.6 and 8.7, the computed value of Xbase ~ Xexit matched the measured 
value within a specified tolerance. 

For a hot-wire measurement point M, (Xm» Ym)> in the boundary layer referenced 
to point P, (Xp, Yp), on the model surface the surface-normal distance, h, is evaluated by 

h = V(Xp-Xm) 2 + {Yp-YmJ 2 . (8.8) 

The appropriate location of P is determined by iteration such that the local slope of 
the model at angle of attack was normal to the slope of the line between points P and M. 
This criterion can be written as 



Since the hot-wire measurements were conducted in a horizontal plane, the angle of 
attack may also be referred to as the sideslip angle p. 
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Fig. 8.7. 






Fig. 8.8. Nozzle and model coordinate systems. 
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